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ABSTRACT 


Gaussian  functions  were  used  as  a  basis  for  the  ex¬ 
pansion  of  the  orbitals  and  correlation  function  for  the 
helium  atom.  The  energy  obtained  with  a  wave  function  of 
orbitals  only  approached  within  0.03%  of  the  Hartree-Fock 
limit.  Except  for  a  slight  rounding  at  the  nuclei  the 
density  distribution  was  identical  to  that  given  by  an  ex¬ 
pansion  in  STF.  Introduction  of  the  correlation  function 
improved  the  energy  markedly  to  within  0.4%  of  Pekeris ' 
value.  The  correlation  function  found  indicated  a  com¬ 
paction  of  electron  density  and  increased  electron  repulsion. 
A  lower  total  energy  was,  however,  achieved  by  an  increase 
in  the  electron  density  at  the  nucleus.  These  results 
were  achieved  in  all  cases  with  a  quite  reasonable  (at 
most  six)  number  of  terms  in  the  expansions. 

The  formalism  of  a  rather  general  SCF  theory  was  ex¬ 
tended  to  cover  more  cases.  Virtual  orbitals  and  ionization 
energies  were  investigated  and  a  more  general,  where 
possible,  formulation  in  terms  of  the  general  SCF  theory 
given.  The  constancy  of  orbitals  between  neighbouring 
states,  which  is  basic  to  certain  uses  of  orbital  energies, 
(i. e. , ionization  energy,  excitation  energy)  was  investi¬ 
gated  for  small  atoms  and,  within  the  Pariser-Parr  context, 
in  polyatomic  molecules.  The  orbitals  were  found  to  be 
similar  in  ground  and  ionized  states  and  the  constancy  be¬ 
came  better  as  the  system  became  larger. 


- 


The  implications  of  the  use  of  the  perturbation 
method  in  describing  NMR  and  EPR  phenomena  were  investiga- 
ted.  The  treatment  was  found  to  lead  to  the  acceptance  of 
concepts,  such  as  coupling  constants  and  well  defined  spin 
states  of  electrons  and  nuclei,  which  obscure  the  funda¬ 
mental  nature  of  the  phenomena.  An  alternative  and 
fundamentally  different  variational  approach  was  developed. 
With  an  illustrative  system  of  two  electrons  and  two  pro¬ 
tons  the  salient  factors  involved  in  the  phenomena  were 
investigated  quite  readily.  A  simple  numerical  example 
indicating  the  relation  of  electronegativity  serves  to 
illustrate  the  ease  with  which  other  relationships  might  be 


demonstrated 
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PREFACE 


While  this  research  has  been  concerned  with  es¬ 
sentially  distinct  problems,  each  with  its  particular 
interest  and  pertinent  results,  the  approach  in  all  cases 
is  pervaded  with  a  desire  to  resolve  problems  for  which 
there  is  no  exact  solution,  with  a  description  which  can 
be  brought  closer  to  the  exact  description  with  the  aid 
of  the  variational  technique  on  some  basis  set. 

In  part  of  this  work  the  representation  of  one- 
electron  orbitals  and  of  two-electron  correlation  functions 
with  Gaussian  basis  sets  is  investigated.  In  another 
section  variational  calculations  for  open-shell  configura¬ 
tions  are  investigated  and  in  a  final  section  in  an 
attempt  to  provide  a  more  satisfactory  description  of 
nuclear  magnetic  resonance  and  electron  paramagnetic 
resonance  phenomena  the  variational  calculation  of  the  spin 
functions  of  electrons  and  nuclei  with  spinors  as  basis 
functions  is  also  considered. 
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PART  ONE 


THE  USE  OF  GAUSSIAN  FUNCTIONS  IN  ATOMIC  WAVE  FUNCTIONS 


CHAPTER  I 


DEVELOPMENT 


3 


1-1  Introduction 

Within  the  orbital  approximation  the  application  of 
the  principle  that  any  well  behaved  function  can  be  repre¬ 
sented  in  terms  of  a  complete  set  of  functions'*"  and  the  use 

2 

of  the  variational  theorem  lead  to  wave  functions  which 
would  more  closely  represent  the  exact  description  of  a 
system.  However,  only  with  the  advent  of  electronic  com¬ 
puters  has  expansion  of  the  basis  set  sufficient  to  yield 
accurate  results  become  feasible. 

Further  improvement  upon  such  an  orbital  description 
of  a  many-body  system  can  be  achieved  by  introducing  inter- 
particle  coordinates  into  the  argument  of  the  wave  function 
through  a  "correlation  function".  Here  again  an  adequate 
representation  of  the  correlation  function  can  be  obtained 
with  an  appropriate  basis  set. 

Clearly  it  is  possible  to  obtain  rather  accurate  des¬ 
criptions  of  systems  given  a  sufficiently  large  basis  set 

and  time  to  perform  the  variational  calculation;  perhaps 

5 

the  best  support  for  this  approach  is  that  Pekeris '  cal¬ 
culated  value  for  the  energy  of  helium  is  now  quoted 
instead  of  the  experimental.  However,  except  for  the  simp¬ 
lest  systems,  mathematical  manipulations  leading  to  the 
necessary  integrals  and  the  evaluation  of  these  integrals 

may  rapidly  become  intractible , *  even  with  the  aid  of  a 

* 

Compare  for  example  the  formulation  of  the  integrals  of 
Roothaan  and  Weiss^  (for  e-mr  basis  functions)  for  helium 
to  those  of  this  work.  The  disparity  is  accented  in  more 
complex  systems. 
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computer,  unless  care  is  exercised  in  the  choice  of  the 
form  of  the  basis  functions. 

The  use  of  Gaussian  functions  as  a  basis  for  the 
variational  calculation  of  orbitals  comprising  wave  func¬ 
tions  of  atomic  and  molecular  systems  has  been  discussed 

S  6 

by  several  authors^’  .  More  recently  the  advantage  of 
representing  two-electron  functions  by  expansions  in 
Gaussian  functions  has  been  noted''7.  In  both  instances  the 
mathematical  manipulations,  although  complex,  lead  to  ana¬ 
lytic  integrals  (up  to  the  error  integral)  and  consequently 
facilitate  the  computational  steps,  thus  admitting  the 
possibility  of  the  solution  of  the  large  problems  of  inter¬ 
est  to  the  chemist. 

Recognizing  this  advantage  of  the  form  of  Gaussian 
basis  functions,  an  investigation  has  been  undertaken  to 
determine  whether  a  reasonable  number  of  Gaussian  functions 
can  yield  wave  functions  with  an  energy  value  comparable  to 
that  obtained  from  other  sources,  and  whether  a  correlation 
function,  which  is  an  expansion  in  a  set  of  Gaussian  func¬ 
tions  having  the  interelectron  distance  as  the  argument, 
can  yield  accurate  energies. 

The  helium  atom  was  used  as  the  test  case.  The  number 
of  basis  functions  used  was  to  some  extent  limited  by  the 
capabilities  of  the  computer  available. 

1-2  The  Wave  Functions 


The  total  wave  function  for  the  non-correlated  case  is 
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given  as 

P  =  A?  =  A  (0  (1)  0  (2)  ) 

where  A  is  an  antisymmetrizer  and  0(1)  and  0(2)  are  spin- 

orbitals  associated  with  spin  and  <^>  respectively.  With 

0  written  as  ^(1)<*(1)  the  orbital  (1)  is  represented 

by  an  expansion  in  Gaussians  as 

^(1)  =  Zp  c.G. (1)  ,  G. (1)  =  exp  (-m.r?)  ; 

i=l  1  1  1  11 

£  is  the  number  of  basis  functions,  and  r^  is  distance  of 

electron  1  from  the  nucleus. 

For  the  correlated  case  the  total  wave  function  is 

given  as 

*  =  A  r.  A  (0(1)  0(2)  %(12)  ) 
where  the  correlation  function  '%(12)  is  expanded  as 

7(12)=Sq  dsGs(12)  ,  G  (12)  «  exp  (“"igifp  ; 
s  =1 

q  is  the  number  of  basis  functions  and  r^  is  the  distance 
between  electrons  1  and  2. 

Since  there  are  only  two  electrons  the  result  of  the 
antisymmetrizer  is  the  same  in  both  cases  and  can  be  writ¬ 
ten  as  the  Slater  determinant; 


A^  = 

0(1) 

0(1) 

0(1) 

0(2) 

^  (12) 

> 

II 

0(2) 

0(2) 

0(2) 

0(1) 

Frequent  reference  will  be  made  to  the  work  of  Roothaan 

4 

and  Weiss7  Their  total  wave  functions  are  formulated  as 
above  except  that  the  expansions  are: 

A»  =  £m  a.p^r)  ,  p,  (r)  =  (  t  r)1  e"  *  r 
i=0  1  1  1 


bdi  ■  i 
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X(t12) 


,n 


i:0  w(ri2) 


yri2} 


There  can,  of  course,  be  no  orthonormality  constraint  on 
the  ^  and  ?C  . 


I~3  The  Variational  Equation 

The  energy  can  be  written  in  bra-ket  notation  as 
E  /H/  '&>/<£'/ 

where  H  =  -E2  (  V^/z  +  —  )  +  —  with  V  ^  denoting 

i=l  1  ri  r12  1 

the  Laplacian  for  the  kinetic  energy.  An  arbitrary  varia¬ 
tion  in  the  energy  is  given  by 

£e  =  [  £  «  jp /H/ jP  ^>  )  -  E  /  <*/§} 

For  the  condition  of  minimum  energy  S E  =  0,  and  since 
<& /&}  must  be  finite  one  obtains 

-  E  =  0 

giving  two  equations  which  are  Hermitian  conjugates 

<3P/K/S&>  -  E  <&/?&>  =  0 
<£3P/k/&>  -  =  o 

The  scalar  products  <iz'/H / £ !py  and  may 

be  simplified.  I  =<'& /K/£JPJ>  =  <  A  V'  /H/  £  A  and  since  the 

variations  are  independent  of  orbital  occupancy  one  has 
I  =<A^/H/AS’^>  =^A^/HA /Spy  .  For  ^  real  and  HA  real 
and  Hermitian  I  =<£  ^/HA/A  ^  >  =  <X^/H/A2  ^>.  But  A2  ^  =KA  p 
where  K  is  a  constant®  so  that  I  =  K  <1T^/H/A^)>  .  Simi¬ 
larly  =  K<^/A  fsy. 
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Then  <  (p  / R/£&>-  K[<T^/H/A^>-  E  ^<f>/H/A  £">  ]  =0 

Since  K/0  this  gives  as  the  variational  equation 

</<T^/H/A^>  -  E <fP/k  =  0 


1-4  The  Pseudo-Eigenvalue  Equations 


The  variation  <7/^/  may  be  writ  ter  for  the  non- 
correlated  case  as  and  for  the  corre¬ 


lated  case  as  </*^/  <p  .  Here 

represents  the  variation  in  the  bra  <  ^ /  when  is  given 
the  arbitrary  variation 

For  the  correlated  case  the  variations  £  and  £X  are 
independent  and  the  variational  equation  can  be  written  as 
two  coupled  equations. 

S/[<ryy  /H/A^>-  E  ~  /A^}l  -  0  1-1 


1-2 


^C<-f^/H'/A^>  "  E<-|^/A^>  3  -  0 

By  virtue  of  the  linear  expansion  form  chosen  for  9^ 
and  ^each  of  the  above  equations  may  be  broken  up  into  a 
set  of  equations  which  in  turn  can  be  represented  by  a 
single  matrix  equation. 

In  the  first  equation  £X<£ /  can  be  written  as 

c)/C 


f  XcJiX-,  /  =  (  §  d  )  \  j and  since  the  £d  are  inde- 

by-  s=i  s  '  s 

pendent,  the  equation  becomes 

ScL  ^<-4rx-  /H/A^>  -  E  4-— -  /A^>  ]  -  o,  S  -  l,q. 


a  d 


d  ft 


Further,  since  each  S  d  is  arbitrary  then  in  each  of  the 

s 

above  equations  the  quantity  within  the  square  bracket 
must  equal  zero. 
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The  scalar  products  can  be  expanded  as 


>Sd  exp[-n  ry  0  0  2  - 

< - 2 - /H/£|,ds,  exp[-ns,r£2]  A  0  0> 

s 

which  on  differentiation  becomes  the  sum  of  scalar  products 
S^,<expC-ngr22]  0  0  /H/  dg,  exp[-ng  .r^]  A  0  0>s  Sd,Hds,  . 

The  two  terms  in  each  of  the  above  equations  become  a  sum 
of  terms  and  one  obtains  the  set  of  equations 

Ss-(Hss'ds-  -  ESss'ds'>  -  0  -  s  ■  1.4, 


or  the  equivalent  form,  the  matrix  equation 


Hd  d  =  E  3d  d 


1-3 


cL  ci 

Ha  and  Sa  are  matrices  with  elements: 


where  d  is  a  vector. 

Hss'=<exp["n^P]  ^1)V^(2)/H/  exp[-nQ,r^p]  ^(1)  y^(2)> 


L  -L  . 

s  12 


12- 


Sss,=<expC-nsri2]  ^  ^  2)/exp[-ng  ,r^2]  ^(l)  j^(2))where  the 
spin  has  been  integrated  out  for  a  spin-independent  hamil- 
tonian 

When  the  variation  <C'|  is  written  as 


a 


a  ^ 


the  matrix  formulation  of  equation  (1-2)  can  be  developed  in 


a  manner  analogous  to  the  above  to  give 

Hc  o'  =  E  Scc* 


1-4 


v  0  0 

where  c  is  a  vector.  H  and  S  are  matrices  with  elements, 
after  integration  of  spin,  as 


9 


,  =^exp[-mir^]  V^2)  X  ( 12)/H/exp[-roi  ,  12)> 

S?it  =<exp[-mir^]/?(2)  X  ( 12)/exp[-mi , r^] ¥tZ)Xl  12) 

Equation(I-4)is  the  matrix  eigen-function  formulation 
of  the  variational  problem  for  the  non-correlated  case. 

The  two  matrix  equations 

Hdd  =  E$1  and  Hcc  =  E§c 

with  the  common  eigenvalue  E,  are  taken  as  the  statements 
of  the  variational  problem  for  the  correlated  case. 


I~3  Formulation  of  the  Matrix  Elements 

On  substitution  for  and  X  the  scalar  product 
matrix  elements  of  the  foregoing  may  be  expanded  and  writ¬ 
ten  in  integral  form.  For  the  correlated  case: 

Hd  1=2.2.(2.2.(c.c.,c.c.,H..,  .... 
ss'  1  i*  J  O'  1  1  J  0  n'jj'ss' 


HT.  ,=  2. 2.  ,2  2  1  c  .  c  .  1  d  d  .  H.  .  ,  .  . ,  , 

11 '  0  0  s  s  0  0  s  s '  n'  00  ss 


S  2.2. ,  2  .2  .  ,  c  .  c  .  , c  .  c  . , S . . ,  .  . .  . 

ss'  1  i'  j  0  1  1  0  a  11  10  ss 


St  .  ,  =  2. 2.,  2  2  ,c  .c.  tdd_  ,S.  , 

11 '  j  j '  s  s '  j  j  s  s  ii'oj'ss' 


where 

H 


ii'jj'ss'"  fexp  [  rairl  mjr2  nsr12] 


2  2  2 

Hexp  [-mi,r1  -  m.(r?  -  ns,r12]  d£ 


Sii' JJ ' ss*  =  lexp  [_mirl  '  mjrI  ■  nsr12] 


2  2  2 

exp  [-m^r^  -  m  ^  , r^  -  ns,r12]  dC 
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and  d  represents  the  volume  element  for  all  space  for  the 
two  electrons. 

For  the  non-correlated  case  there  is  only  the  one 

C-*  c  — * 

matrix  equation,  He  =  ES  c ,  and  the  elements  reduce  to 

H  •  ■  •  —  2  « IS  •  •  c  •  c  •  *  H  •  *  I  •  *  *  .  3  •  ■  §  “ *  S  •  S  •  •  c  •  c  •  •  3  •  <  ■  «  •  • 
ii  J  J  J  ,  1  n'n'  ’  ii '  ; i  '  n'n.i' 


where 


a  2  2  o 

Hii'aj,=  )exP["miri  "  mjr2-*  H  expC-nn.r^  -  m^r^]  d  u 

r  2  2  2  2 
Sii'jj,=  )exPt“mj_ri  ~  mjr2-*  exp  ”  mj  i r2^  d  ^  * 


1-6  Formulation  of  the  Integrals 


The  matrix  elements  of  the  foregoing  lead  to  integrals 
of  the  form 

I  =  [exp  [-Sja-jT?  -  jbijrij ] 

Mexp  C-2kakrk  -  skibkirki^ 

where  M  is  some  operator. 

6  7 

Boys^* '  has  developed,  in  a  series  of  papers,  a  gen¬ 
eral  method  of  evaluating  integrals  of  this  form.  Although 
much  of  this  method  depends  upon  the  execution  of  the 
proper  mathematical  manipulations  at  the  proper  time,  there 
is  in  addition  the  introduction  of  the  density  kernel,  which 
allows  both  the  very  elegant  formulation  of  the  integrals, 
and  subsequently  the  ease  of  their  solution. 

The  density  kernel  allows  an  integral  <  ®1/Q/*2>,  where 
$  is  a  function  of  many  electron  co-ordinates  and  ^  is  a 
sum  of  one  electron  operators,  to  be  expressed  as  a  sum  of 
integrals  ^  where  the  0^  are  functions  only  of 
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the  electron  involved  in  Q^.  Thus  the  density  kernel 
operator  for  electron  i  involves  integration  over  the  co¬ 
ordinates  of  all  other  electrons  and  retains  the  identity 
of  the  coordinate  of  electron  i  as  belonging  either  to  the 
bra  or  to  the  ket.  This  last  requirement  is  necessary  since 
certain  of  the  operators  will  be  differential  operators 
and  act  only  on  the  coordinates  of  r^  in  either  the  bra  or 
ket  but  not  both. 

Letting  P(x'/x" )  denote  the  operator  replacing  x"  by 
x',  then  the  density  kernel  operator  R(r'>~r,,,i)  for  electron 
i  is  defined  by 

^0-i /R(r '  ,r" , i)/0?>  =  (tcN  $  dr.)  P(r '/r.  )0,  P(r"/r.)0? 

1  2  d/i  J  11  12 

=  c^Cr'  r") 

so  that  this  operation  denotes  coordinates  of  electron  i 
coming  from  the  bra  as  r'  and  those  from  the  ket  as  r"  and 
integrates  over  the  others  to  give  the  factor  c^.  Then 
the  integral 

<01/SiQi/02^  =  ^drP ( r/r '  )  P(r/r")Q(r,,)Si  <  01/R(r;r" ,i)/02>  . 

It  will  be  noticed  that  the  operator  has  been  moved 
to  the  left  of  the  coordinates  of  the  bra  and  ket  but  since 
the  density  kernel  operator  has  denoted  the  coordinates  of 
electron  i  coming  from  the  bra  as  r'and  those  coming  from 
the  ket  as  r"  then  the  operator  need  only  be  expressed 
as  Q(r")  in  order  that  it  operate  correctly.  Once  the 
operation  by  Q  has  been  performed  there  is  no  longer  any 
need  for  distinction  between  r'  and  r"  and  so  the  operators 
P(r/r ' )  and  P(r/r")  operate  on  r'  and  r"  leaving  the  integral 
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JdxF(*)  to  be  evaluated.  F(r)  is  the  result  of  Q(r") 

operating  on  ^(r^ r" )  followed  by  operation  by  the 

P's,  i. e. , F  (r)  =  P(r/r '  )  P(r/r")  Q(r")  Sicif  i(r>'r")  . 

The  utility  of  this  form  can  be  appreciated  when  it  is 

2  1 

noted  that  the  integrals  for  several  operators  Q.=-V.,  — 

l  1  ri 

etc.,  can  be  evaluated  using  the  same  density  kernel. 

7 

Boys"  general  approach  is  particularized  to  the  helium 
system  in  Appendix  I.  When  the  hamiltonian  is  separated 
into  kinetic  energy,  nuclear  attraction  and  repulsion  of  the 
two  electrons  the  integrals,  as  developed  in  Appendix  I, 
have  the  forms  listed  below.  It  is  worthwhile  emphasizing 
that  the  integrals  listed  are  all  analytic. 


H, 


ii' 11  ss 


, =  H(KE1)  +  H(KE2)  +  H(NA1)  +  H(NA2) 


ii' JO 'ss' 

+  H(Rep) 
ii' jj 'ss' 


ii ' j j ' ss 


ii ' j  j ' ss ' 


ii ' JJ ' ss ' 


H(KE1) 
ii' j j ' ss 


_  „  5(albl  +  bl°l  +  °lal)  3/2 
■ 


H(KS2)  =  Gp 
ii  '  j  j  '  ss  '  (a2  +  b2) 


^  a2^2+  ^2c2+  c2a2^  3/2 

- 572 -  71 


H(NA1) 

11  JJ  SS' 


2mG. 


ai+  bi 


2tcG  p 

ii  j  j  ss  2  < 


H(Rep) 
ii' j j ' ss ' 


2tu5//2(U  +  V)’1/2 

uv  +  uw  +  vw 


Sii'ad'ss'=  (uv  +  +  )'3/2  *5 
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where : 

G-,  =  T[^//^(m.+  ra  . ,  +  n  +  n  ,) 

1  j  j  s  s 

Gp=  n^//^(m.+  m.  ,+  n  +  n  ,) 

<—  11  s  s 

n  (n  +  n  , ) 
s  s  s 

a,  =  m.  +  n  -  —  ■  .  — y-- ■— t . _ — 

1  l  s  m  .  +  m  . ,  +  n  +  n  , 

3  0  s  s' 

n  (n  +  n  ,  ) 
s  s  s 

a0=  ra.+  n  -  — - - — - r - 

2  j  s  ra .  +  m .  ,  +  n  +  n  , 

u  i  i '  s  s ' 

n  ,  (n  +  n  ,  ) 
s '  s  s '  ' 

m .  .  +n  .  — — - — - — - 

l  ’  s'  m  .  +  m  . ,  +  n  +  n  , 

3  0  s  s 

n  , (n  +  n  , ) 

_  s'  s  s 

m  j  '  n  s  '  m .  +m .  ,  +  n  +  n  , 

°  1  l  s  s 

n  n  , 

=  s  s' _ 

1  ra  .  +  m  . .  +  n  +  n  , 

3  o'  s  s' 

n  n  . 

c  =  s  s' _ 

2  mi+  nr  , +  ng+  ns , 

U  =  m .  +  m .  ,  V  =  ra  .  +  m  .  ,  W  =  n  +  n  , 
li  GO  s  s 

I~7  Self-Consistent  Field  Solutions 

A  basis  program  to  calculate  the  matrix  elements  was 
written  for  the  IBM  1620  and  7090. 

Solutions  of  the  matrix  equations  for  the  eigenvectors 
and  eigenvalues  were  carried  out  by  performing  a  Schmidt 
Orthonorraalization  followed  by  a  Jacobi  diagonalization  and 
a  search  for  the  lowest  eigenvalue  and  eigenvector.  These 
steps  are  examined  in  detail  in  Appendix  II. 


Using  the  basis  program  in  conjunction  with  programs 
for  the  procedures  mentioned  above,  a  master  program  was 
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constructed  which  solved  the  matrix  equations  given  by  a 
set  of  basis  functions  and  trial  vectors  for  a  new  set  of 
vectors.  These  last  vectors  then  became  the  trial  vectors 
and  the  calculation  of  the  matrix  elements  and  solution  of 
the  matrix  equations  was  repeated  until  the  output  and  in¬ 
put  vectors  were  the  same  to  five  figures.  The  vectors  and 
energy  found  at  this  point  were  taken  to  be  the  self- 
consistent  field  (SCF)  values  for  a  given  basis  set.  Figure 
1  presents  a  flow  diagram  of  the  master  program  with  associa¬ 
ted  programs. 

1-8  Variation  of  the  Basis  Set 

Variation  of  the  exponential  parameters  m  and  n  of  the 
Gaussian  basis  sets  was  attempted  by  two  associated  pro¬ 
grams.  The  first  of  these,  referred  to  as  "Mavle",  varied 
one  parameter  at  a  time  until  a  temporary  best  value  of 
that  parameter  was  found  and  then  continued  on  to  the  next 
parameter,  and  so  on  returning  finally  to  the  first  para¬ 
meter  and  repeating  the  cycle  until  the  best  value  of  all 
parameters  had  been  found.  The  second,  the  method  of  quad¬ 
ratic  minimization, is  a  method  of  solving  for  the  best 
exponential  parameters  when  the  energy  surface  in  the  space 
of  the  exponential  parameters  had  been  mapped  in  the  region 
of  the  minimum.  The  first  procedure  is  examined  in  detail 

in  Appendix  III  and  the  second  has  been  discussed  in  some 

g 

detail  by  Ransil.  Of  the  two  methods  the  former,  somewhat 
less  elegant  procedure,  was  found  to  be  the  most  satisfac¬ 
tory.  Additional  provision  was  made  for  varying  the  basis 
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set  by  changing  the  number  of  terms  in  either  or  both  of 
the  orbital  and  correlation  expansions.  The  storage 
capacity  of  the  machine  limited  this  flexibility. 

I~9  Electron  Density 

Three  forms  of  the  electron  density  were  used  in  an 
effort  to  provide  insight  into  the  detailed  character  of 
the  wave  function. 

Defining  the  Radial  Electron  Density  as 

RD(r2)  =  [  ^2d  t-L  sin  62d©2d02]  r| 

gives  RD(r2)  =  G24tc  r2  where  G2  is  the  density  kernel  for 
electron  two. 

The  Point  Electron  Density  is  defined  as 
PD(r2)  -  -  G2 

The  appropriate  form  of  G2  ,as  given  in  section  AI-1, 
was  inserted  into  these  formulae  and  a  program  written  to 
calculate  the  densities  for  various  values  of  r. 

The  program  could  be  modified  to  provide  ,r2)/2 

allowing  the  probability  function  to  be  investigated  for 
particular  values  of  both  electron  coordinates. 
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Figure  1 

Flow  Diagram  of  Master  Program 


End 


CHAPTER  II 


RESULTS  AND  DISCUSSION 
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II-l  The  Orbital  Representation 

Table  I  displays  the  best  exponential  parameters, 

SCF  coefficients,  and  energies  found  for  the  non-correlated 
wave  functions  when  the  orbitals  are  given  by  one~to  six- 
term  expansions  in  Gaussians. 

A  comparison  of  the  energies  given  in  column  nine 
indicates  that  further  expansion  of  the  basis  set  would 
improve  the  energy  only  slightly. 

The  best  energy  achieved  here  differs  by  only  0.03% 
from  the  Hartree-Fock  energy  of  -2.862  au. *,  which  is 
in  addition, the  best  energy  of  Roothaan  and  Weiss  using  a 
five-term  expansion  in  restricted  Slater  functions. 

In  view  of  this  one  may  conclude  that  a  reasonable 
number  of  terms  (five  or  six)  in  a  Gaussian  expansion  gives 
energies  for  the  helium  atom  which  are  close  to  those  using 
Slater  functions  and  which  compare  favourably  to  the  Hartree- 
Fock  energy. 

The  radial  density  distributions  and  the  point  density 
distributions  are  displayed  in  Tables  II  and  III  and  figures 
2  and  3  for  several  of  these  expansions  and  for  the  single- 
and  four-term  Slater  expansion  of  Roothaan  and  Weiss^ 

It  is  apparent  that  the  effect  of  increasing  the  number 

of  terras  in  the  Gaussian  expansion  from  three  to  six  is,  as 

viewed  from  the  point  density  plot,  to  raise  the  electron 

^Throughout  this  thesis  au.  refers  to  Hartree  Atomic  Units 
(i.e.  1  au.  =  27.204-  ev) . 
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density  in  the  region  of  the  nucleus  (with  a  slight 
increase  in  the  slope),  and,  from  the  radial  distribution 
plot,  to  move  the  maximum  somewhat  closer  to  the  nucleus  - 
clearly  an  attempt  to  correct  for  the  zero  slope  of  the 
Gaussian  form  at  the  nucleus.  On  the  other  hand  the  in¬ 
crease  in  the  number  of  terms  in  the  Slater  expansion  is 
effecting  an  increase  in  the  density  about  the  nucleus 
(with  slight  decreases  in  the  slope)  and  flattening  the 
radial  distribution  curve  -  in  this  case  an  attempt  to 
correct  for  rather  peaked  character  of  the  Slater  function. 
However,  the  distributions  as  given  by  the  best  Slater  and 
best  Gaussian  expansions  are  essentially  the  same,  an  in¬ 
dication  that  although  the  Gaussian  expansion  is  approaching 
the  representation  of  the  orbital  from  a  different  direction 
than  the  Slater  expansion,  the  result  is  the  same. 

One  might  have  expected  the  inadequacy  of  the  Gaussian 
form  in  the  region  of  the  nucleus  to  preclude  a  good  IS 
orbital  representation  or  indeed  even  an  approach  to  the 
Slater  representation.  However,  in  view  of  these  results, 
one  may  suggest  that  a  Gaussian  expansion  can  represent 
adequately  the  IS  type  of  orbital  with  a  reasonable  number 
of  terms. 

In  regions  further  from  the  nucleus  the  Gaussian  form 

would  be  more  plausible^  and  representations  of  2S,  2P- 

2  2 

type  functions  with  basis  functions  of  the  form  r  exp[-mr  ] 

p 

and  z  exp[-nr  ],etc.,  should  be  practicable  -  certainly  as 
feasible  as  with  Slaters  -  so  that  Gaussians  give  promise 
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of  an  adequate  representation  of  the  orbitals  for  many- 
electron  problems. 

II-2  The  Correlation  Function  Representation 

Correlation  functions  given  by  two-,  four-,  and  six- 
term  expansions  were  associated  with  one-electron  protions 
generated  by  one-,  three-,  and  six-term  expansions  to  give 
the  correlated  wave  function.  Table  IV  displays  the  best 
exponential  parameters,  SCF  coefficients  and  the  energies. 

The  best  energy  obtained  is  0.4%  above  Pekeris'^ 
value  and  0.3%  above  the  best  value  of  Roothaan  and  Weiss 
with  their  correlation  function  for  a  closed  shell  function. 

A  comparison  of  the  energies  for  p=  2;  q=  2,4,  and  6  in 
the  table  indicates  that  extension  of  the  basis  set  for  the 
correlation  function  beyond  the  first  two  terms  improves 
the  energy  only  very  slightly. 

It  is  clear  then,  that  an  improvement  in  the  energy 
(some  71%  of  the  difference  between  Hartree-Fock  and  experi¬ 
mental),  comparable  to  that  using  other  functions  can  be 
produced  by  a  correlation  function  given  by  a  quite  small 
expansion  in  Gaussians. 

In  all  cases  the  introduction  of  the  two-electron 
variable  results  in  an  improvement  in  the  energy.  The  source 
of  this  improvement  may  now  be  investigated. 

The  function  given  in  the  tables  is  depicted  sche¬ 
matically  in  figure  4.  The  classical  form  of  the 

12 

function  of  Hartree  and  Ingman  which  has  a  minimum  for 
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r-^2=0  and.  a  finite  limit  for  ^-^2*  °°? and-  "X  function 
determined  by  Roothaan  and  Weiss  are  given  for  comparison 
in  the  same  figure. 

The  Gaussian  correlation  function  does  not  introduce 
a  lowering  of  weight  of  the  wave  function  on  close  approach 
of  the  two  electrons  which  both  the  Roothaan  and  Weiss  and 
the  Hartree  and  Ingman  functions  would  accomplish.  It 
does,  however,  approach  an  asymptotic  value  as  r-^-*  00  > 

as  required  classically,  while  the  Roothaan  and  Weiss  func¬ 
tion  goes  to  ~  00  0 

Although  both  the  Roothaan  and  Weiss  and  the  Gaussian 
correlation  functions  imply  a  quite  non-classical  behaviour, 
the  inadequacy  of  the  former  would  not  necessarily  be  impor- 
tarfc  for  the  over  all  description  since  the  orbital  parts 
would  certainly  cause  the  total  function  to  go  to  zero  as 
r12*  00  suggesting  a  well  behaved  total  wave  function.  Such 
a  happy  resolution  of  the  problem  is  not  possible  in  the 
Gaussian  case. 

The  behaviour  of  the  total  wave  function  with  Gaussian 
functions  can  be  investigated  by  a  comparison  of  the  curves 
for  the  correlated  and  non-correlated  functions,  for  the 
variation  of  /i^(r-^,  r^=  0.6)/  as  the  coordinate  r^  is 
changed  along  a  line  through  electron  two  and  the  nucleus. 

As  figure  5  shows  there  is  an  increased  probability  for 
the  first  electron  in  the  region  of  the  second  compared  to 
the  probability  for  the  same  radial  distance  on  the  oppo¬ 
site  side  of  the  nucleus.  Indeed  the  most  probable 


. 

nr  rvnl 
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position  of  the  first  electron  is  between  the  nucleus  and 
the  other  electron,  not  as  one  might  have  expected,  on  the 
opposite  side.  Clearly  the  introduction  of  explicit  depen¬ 
dence  on  the  interelectron  distance  with  Gaussians  is 
inducing  a  compactness  of  the  electrons  -  a  behaviour  which 
is  imposed  essentially  by  the  ^-function  alone. 

Such  a  compaction  and  accompanying  increased  repulsion 
can  result  in  an  overall  lowering  of  the  energy  if  the  build¬ 
up  can  be  made  to  occur  closer  to  the  nucleus.  Comparison 
of  the  energies  between  Tables  I  and  IV  shows  that  the  improve¬ 
ment  in  the  energy  with  correlation  is  most  marked  the  larger 
the  expansion  for  the  one-electron  part.  But,  as  the  curves 
of  figures  3  and  6  show,  the  larger  the  expansion  the  greater 
the  density  at  the  nucleus,  hence  a  greater  improvement  on 
correlation.  Further,  comparison  of  the  exponents  of  the 
one-electron  parts  of  the  correlated  and  non-correlated 
functions  as  given  by  Tables  I  and  IV  indicates  a  tendency 
for  an  increase  in  the  electron  density  for  the  correlated 
case  in  the  region  of  the  nucleus  -  a  behaviour,  at  least 
for  the  orbital,  quite  the  opposite  to  that  in  Roothaan  and 
Weiss'  best  function.  It  might  be  pointed  out  here  that 
for  the  poorer  functions  of  Roothaan  and  Weiss  a  similar 
behaviour  (the  orbital  exponent  increase)  occurs  with  cor¬ 
relation  -  an  indication  that  condensation,  even  with  their 
function,  can  be  a  first  step  toward  energy  improvement. 

It  appears  then  that  the  correlation  obtained  here  is 
not  toward  a  decrease  in  the  weighting  of  the  function  for 
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rl2*  0.  Rather  it  is  an  attempt  to  improve  on  the  energy 
by  a  compensatory  process  -  the  compaction  of  the  electrons, 
and  resulting  repulsion,  being  compensated  for  by  a  greater 
nuclear  attraction. 

That  the  improvement  in  the  energy  must  be  produced  by 
such  non-classical  correlation  may  be  due  to  the  Gaussian 
form  of  either  the  one-electron  functions  or  the  correlation 
function  itself.  The  investigation  of  a  Gaussian  correlation 
function  with  one-electron  functions  given  by  Slaters  would 
certainly  shed  some  light  on  this  aspect. 

Although  the  source  of  improvement  in  the  energy  is 
questionable  one  may  suggest  that  improvement  in  the  energy 
with  a  correlation  function  given  by  an  expansion  in 
Gaussians  is  feasible. 
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Table  I 


Parameters*  for  the  expansion  for  the  non-correlated  case. 
Parameter  Parameter  position  in  the  expansion  Energy 


type 

1 

2 

3 

4 

5 

6  (au.) 

1 

c 

0.263 

m 

0.76  7 

-2.301 

2 

c 

0.366 

0.581 

m 

0.532 

0.410 

-2.747 

3 

c 

0.228 

0.490 

0.406 

m 

0.383 

2.00 

13.6 

-2.836 

4 

c 

0.143 

0.389 

0.413 

0.265 

m 

0.293 

1.21 

5.59 

37.0 

“2.855 

5 

c 

0.090 

0.302 

0.386 

0.304 

0.177 

m 

0.235 

0.820 

3.06 

13.4- 

88.6 

-2.860 

6 

c 

0.052 

0.215 

0.335 

0.326 

0.226 

0.128 

m 

0.189 

0.564 

1.77 

6.14 

26.3 

170.0  ”2.861 

*  The  energy  to  four  figures  is  insensitive  to  changes  of 
two  to  seven  parts  in  the  third  figure  of  the  exponential 
parameters  m  and  to  one  part  in  the  third  figure  in  the 
linear  coefficients  c. 


,  . 

Q.<?€ 
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Table  II 


Radial  distribution  of  electron  density 
for  the  non-correlated  case. 


Radius 

(a0) 

One 

Optimized 

Slater 

Roothaan 

and 

Weiss 

Gaussian 

3 

Expansions 

6 

0 

0 

0 

0 

0 

0.1 

0.137 

0.152 

0.142 

.152 

0.2 

- 

0.413 

0.418 

0.414 

0.3 

0.628 

0.639 

0.635 

0.637 

0.4 

— 

0.786 

0.764 

0.783 

0.5 

0.889 

0.856 

0.844 

0.857 

0.6 

0.914 

0.865 

0.877 

0.866 

0.7 

0.888 

0.832 

0.860 

0.831 

0.8 

- 

0.772 

0.800 

0.772 

0.9 

- 

0.699 

0.712 

0.700 

1.0 

0.658 

0,621 

0.  619 

0.625 

1.1 

— 

0.543 

0.533 

0.547 

1.2 

— 

0.469 

0.457 

0.473 

1.4 

0.326 

0.341 

0.340 

0.345 

1.6 

- 

0.240 

0.252 

0.  244 

1.8 

0.094 

0.166 

0.181 

0.163 
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Table  III 


Point  distribution  of  electron  density 
for  the  non-correlated  case. 


Radius 

(a0) 

One 

Optimized 

Slater 

Roothaan 

and 

Weiss 

Gaussian 

3 

Expansions 

6 

io"4 

— 

1.741 

1. 261 

1.628 

0.05 

— 

1.463 

1.226 

1.477 

0,1 

1.008 

1.  210 

1.126 

1.210 

0.2 

— 

0.834 

0.836 

0.834 

0.3 

0.  558 

0.  562 

0.561 

0.562 

0.4 

— 

0.391 

0.380 

0.389 

0.5 

0.  282 

0.273 

0.  268 

0.273 

0.6 

0.  202 

0.101 

0.194 

0.191 

0.7 

0.144 

0.135 

0.139 

0.135 

0.8 

— 

0.096 

0.099 

0.096 

0.9 

- 

0.069 

0.070 

0.069 

1.0 

0.052 

0. 050 

0.049 

0.050 

1.1 

— 

0.036 

0.036 

0.036 

1.2 

- 

0.026 

0.026 

0.026 

1.4 

0.013 

0.014 

0.014 

0.014 

1.6 

— 

0.00  7 

0.007 

0.007 

1.8 

0.002 

0.001 

0.001 

0.001 
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Table  IV 


Parameters*  for  Gaussian  expansion  for  the  correlated  function. 


Parameter^  Parameter  position  in  the  expansion  Energy 

-\rr  o  -z  It  c.  C. 


p 

4 

type 

1* 

2 

3 

4 

5 

6 

(au. 

1 

c 

G.  334 

m 

0.818 

2 

d 

87.7 

-84.7 

-2.320 

n 

0. 

0.005 

3 

c 

0.082 

0.174 

0.144 

m 

0.367 

2.0 

13.9 

2 

d 

6.  51 

1.91 

-2.855 

n 

0. 

0.35 

3 

c 

0.082 

0.170 

0.144 

m 

0.369 

1.98 

13.5 

4 

d 

6.32 

0.882 

0.725 

0.695 

"2.857 

n 

0. 

0.15 

0.35 

0.55 

3 

c 

0.080 

0. 166 

0.139 

m 

0.371 

2.00 

13.6 

6 

d 

0.  52 

0.839 

1.03 

0.681 

0.304 

-0.753 

-2.859 

n 

0. 

0.83 

0.33 

0.64 

1.2 

6.8 

6 

c 

0.106 

0.199 

0.297 

0.317 

0.  244 

0.144 

m 

0.239 

0.692 

1.79 

5.56 

22.2 

147. 

2 

d 

0.74-6 

0.  278 

-2.891 

n 

0. 

0.34- 

*The  sensitivity  of  the  energy  to  the  parameters  is  the  same 
for  the  c  and  ra  as  in  Table  I,  for  the  d  and  s  the  energy 
is  somewhat  less  sensitive. 

^Although  complete  flexibility  was  allowed  in  the  determina¬ 
tion  of  the  parameter  n-^  the  best  value  in  all  cases  was 
effectively  zero. 
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Table  V 


Probability 

distribution  of  / jp ( 

p 

r1?  r2=  0.6)/  . 

Gaussian 

expansion 

rl(aG} 

Correlated 

P=  3,  q=  2 

Non-correlated 
P  =  5 

-0.7 

0.025 

0.027 

-0.5 

0.051 

0.052 

-0.5 

0.  Ill 

0. 109 

-0.1 

0.254 

0.219 

0.0 

0.267 

0.245 

0.1 

0.261 

0.219 

0.5 

0.125 

0.109 

0.5 

0.060 

0.052 

0.7 

0.052 

0.027 
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Table  VI 


Point  density  distribution  for  the  correlated  function,, 


Radius 

(aQ) 

Gaussian 

p  6 

q  2 

expansions 

3 

2 

0 

1.700 

1.317 

0.1 

1.258 

1.198 

0.2 

0.859 

0,801 

0.3 

0.  584 

0.  580 

0.4 

0.402 

0,390 

0.3 

0.280 

0,274 

0.6 

0.195 

0,197 

0.7 

0.139 

0.141 

0.8 

0.097 

0,099 

1.0 

0.049 

0,049 

1.2 

0.025 

0.025 

1.4 

0.013 

0.013 

1.6 

0.007 

0.007 
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Figure  2 


Radial  distribution  of  electron  density  for 
the  orbital  of  non-correlat ed  functions 

Legend: 

Roothaan  and  Weiss  -i 
Six-term  u-aussian  I 

One  optimized  Slater  _ _ 


Three-term  Gaussian 
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Point  distribution  of  electron  density 
for  the  orbital  of  non-correlated  functions 

Legend : 

Roothaan  and  V/eiss 


Six-term  Gaussian 
One  optimized  Slater 
Three-term  Gaussian 
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Figure  4 


Dependence  of  ^  upon  the  interelectron  distance 


12 


Legend : 

I  Hartree  and  Ingman 
II  Gaussian 


III  Roothaan  and  Weiss 
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Figure  5 


Comparison  of  the  probability  distribution 

/?(r1}0,  6)/2 

for  the  correlated  and  non-correlated  function 


Legend : 


Correlated  p  =  3>  q  =  2 
Non-correlated  p  =  3 
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Comparison  of  point  density  distribution  for  correlated  and 

and  non-correlated  functions 


Legend : 

Correlated  p  =  6,  q  =  2 
Non-correlated  p  =  6 


PART  TWO 


OPEN-SHELL  CONFIGURATIONS 


CHAPTER  III 


THE  GENERAL  SCF  THEORY 
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III-l  Introduction 

The  determination  of  the  electronic  energy  by  a 
variational  calculation  has  been  restricted  to  those 
systems  that  the  available  SCF  formalism  was  capable  of 
handling  -  namely  closed-shell  and  some  open  shell  con¬ 
figurations^*^.  As  a  result  the  energies  of  open-shell 
systems  have  been  approximated  by  calculations  upon  the 
nearest  tractible  form.  For  example^the  ionization 
potential  has  been  approximated  by  the  energy  of  the  high¬ 
est  occupied  orbital  in  the  parent  closed— shell  ground  state; 
even  some  attempts  at  describing  excited  states  have  been 
made  by  the  use  of  the  virtual  (non-occupied  in  the  ground 
state)  orbitals. 

Recently  more  general  methods  have  been  advanced,  in 

17 

particular  Huzinaga's  '  method  in  which  a  SCF  improvement 

of  one  orbital  at  a  time  is  carried  out,  and  that  of  Birss 
1  P 

and  Fraga  which  formulates  the  SCF  problem  for  a  simul¬ 
taneous  variation  of  all  orbitals. 

Much  of  the  next  three  chapters  is  dependent  upon  the 
latter  authors'  SCF  formulation  and  a  resum£  of  their 
approach  with  the  pertinent  equations  follows. 

III-2  The  General  SCF  Equations 

A  separation  on  the  basis  of  symmetry  designation  is 
made  at  the  outset  which  allows  the  energy  to  be  written  as 
a  sum  of  two  parts:  the  first  involving  orbitals  of  a 


.. 
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certain  species  \i  and  subspecies'*;  the  second,  desig¬ 
nated  below  as  C,  containing  all  terras  which  do  not  involve 
these  orbitals.  Thus  one  has  equation  (4)  of  the  reference 

as : 


E  =  2E .  fW  +  E.E  .f^fVl^’v^+  2E.E  S  E.  ,  f^f  VI^  ’  V(S  +  C 
1  ii  l  3  l  3  13  l  m  v  /S l  ra  im 

with  the  i,j  summations  extending  over  all  occupied  orbitals 

of  symmetry  designation  \i ,  *  and  the  m  over  all  occupied 

orbitals  of  representations  covered  by  v,^  .  The  index  ($> 

cannot  represent  ex.  when  v  =  (i  .  The  fractional  occupancy 

is  represented  by  f  with  appropriate  super-and  sub-scripts. 


The  I 


|ix,  v(3 


ira 


represent  combinations  of  two-electron 


integrals  and  are  given  in  the  reference  as 


im  1  ira  im  im  "im 

the  and  Khl<*‘»v^  being  the  usual  coulomb  and  exchange 

im  im  ° 


integrals 


^  <f\  ™Yl  (2)412/^i  (Dfm  (2)> 


KU<,^  ,(5Pji  (1)^v  (2)4  /fl  (2)> 

Variation  in  the  orbitals  leads  to  the  Hartree-Fock 
equations 


F^ 

l 


<x 


(X 


where 

=  fM-H  +  S.  +  E  E  E  /  f^fvIV^ 

l  l  ii  33  mv^x  imra 


III-l 


Thence,  with  the  introduction  of  the  coupling  operator 

which  effectively  removes  the  off-diagonal  Lagrangian 


- 
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multipliers,  the  set  of  pseudo-eigenvalue  equations 
s|i«y>£«  is  obtained. 

R11^  =  SkrkK  •  III- 2 


-  z 


/  U?  \X<K 

v  ri 


+ 


<°V 


><n^/ 


III-3 


III~3  The  LCAO  Represent ation  of  the  SCF  Equations 

19 

In  a  more  recent  work  'by  the  same  authors  the  extension 
has  been  made  to  the  Linear  Combination  of  Atomic  Orbitals 
(LCAO)  context.  The  orbitals  are  expanded  in  terras  of  a 
suitable  set  of  basis  functions  Xa  to  give 

ir 

=  2  7^  C ,J' .  =  where  the  X  ^  is  a  row 

p  pi  '  1 

vector  whose  elements  are  the  basis  functions  of  symmetry- 
designation  \xoc  and  C^  is  a  column  vector  with  elements  cX  . 

The  matrix  representation  of  the  pseudo-eigenvalue 
equations  is  then 

qM-  -  eM-f  qM  III-4 

—  1  11  —  1 

M 

and  the  representation  of  Rr  is  developed  in  the  reference 


as 


Rq*  =  z\lol  +  z\uk  + 
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with  =  U^S'^  -  s^v^si1*  +  -^-S^W^S^ 

where  #*«  S,u£"  ;  tff=  f£*d£ 

—  k— V  ’  — k  — k  — k 


yB*- 


#*  =  SkW^;  W^=  J2v— k* 


k-k 


The  and  ,  are  respectively 5 the  density  matrix 


=  c£c^+ 

— k  k  k 


(which  may  he  expanded  as  the  supervectors 


I)£)  and  the  representation  of  the  operator  F,  equation 
(Ill-l)^in  the  basis  set. 

An  element  of  _F^ is  given  as 


k 


k,  so 


+  s  r  K^p^Byk)]  - 

st  v,(5u^st  &st  ^k 

and  v  ^  are  supermatrices  with  elements  as  de¬ 

fined  in  the  reference  as 

(^}  dva 

(?k)  dvp  dv0 


and 


tv(p) 

-k 


and  B  are  column  supervectors  defined  as 


a^(4)=  22  fva£’v  Dv 
Xk  m  m  km  .~m 


v 


B^(p)  _  2 s  fvb!^>v  D 
k  m  m  km  s^m 


where  the  m  summation  extends  over  all  occupied  orbitals 
of  symmetry  designation  v  $  . 

The  total  energy  can  then  be  developed  as 

e  =  s.  (f^r  ♦  e^)  . 

ip°*  1  1  11 


III-5 


. 


CHAPTER  IV 


FEATURES  OF  THE  GENERAL  SCF  FORMULATION 
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IV-1  Extension  of  the  Formulation 

The  definition  of  the  two  electron  integrals  as 

2a^J^’v^  -  b^vK^*’ was  intended  by  the  authors 

lm  ini  im  1m  im  ° 

to  cover  all  except  2  states.  However,  it  can  be  shown  to 
be  inadequate  for  configurations  such  as  the  singlet 
,  e .  g . , 1  ( 2Px ,  2Py). 

The  contribution  of  the  electronic  repulsion  to  the 
total  energy  can  be  written  for  the  ^(2Px,  2Py)  configura¬ 
tion  as  Eg  =  y  +  jrPxPy.  there  being  no  contribution 

from  the  JPxPx,  KPxPx,  or 

Since  the  JPxPx  ^  jfxPy  ^  c^ear  that  the  coef- 
ficient  required  (here  a  =  C  )  for  the  integral  J  is 


different  than  that  required  for  the  J' 


•PxPy 


The  situa¬ 


tion  is  similar  for  the  K^x^x  and  KPxP^".  Hence  the  a  and 

b  which  were  superscripted  toward  symmetry  species  only, 

must  also  be  superscripted  toward  the  subspecies. 

The  corrected  defintion  of  would  then  read 

lm 

Tp«,v(?  _  ?  jUSW?  - 

im  im  im  im  ira 

With  this  definition  the  derivation  as  given  by  Birss 

and  Eraga  can  proceed  unchanged  to  their  final  formulation 

IV- 2  The  Ionization  Energy 


The  ionization  energy  may  be  defined  as  E  -  E,  where 

S 

E^  and  E^  are  the  total  electronic  energy  of  the  ground  and 
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20 

ionized  states  respectively*  Koopmanb  theorem  approxi¬ 
mates  E  -  E.  by  the  orbital  energy  of  the  electron  removed, 
§  i 

on  the  assumption  that  the  orbitals  are  unchanged  in  the 
ground  and  ionized  states. 

The  expression  for  Koopmante  theorem  for  a  closed 
shell  ground  state  in  the  formalism  of  Birss  and  Braga  may 
be  readily  found  by  properly  relating  some  of  their  formu¬ 
lae  to  earlier  work. 

16 

Roothaan  formulates  the  total  energy  for  a  closed 
shell  as  E(F)  =  E.E  E  (H^*+  G  where  H^(l)  =  </.  ^ 

and  €r^  the  orbital  energy;is  given  by  . 

From  the  equation  (24)  of  reference  (18)  one  has  for 
the  closed  shell  that 


B =  Tpf-V*  =  e^f.^.But  R/,^=  so  that 

ZL»  3.  3.  3.  3. 3>  3> 

€  ^**  =  6^.  for  the  closed  shell. 

i  li 

Equation  (III-5)  for  the  total  energy  then  becomes 
E(R)  =  E.E  E  (©^*  +  )  =  E.E  E  (6^  +  E^  ) 

'  llloC^ll  1  y  1  u  ^  1  1 


One  then  obtains  E  (R)  =  E  (F) ,  -as  one  must;  an  ex¬ 

plicit  statement  (not  given  by  the  authors)  that  the  total 
energies  by  the  two  methods  are  identical  for  the  closed 

shell. 

When  the  orbitals  pf  the  ground  and  ionized  states 

are  the  same , Roothaan'*' ^  has  shown  that  E  (F)  -  6  =  E .(F). 

g  v  v 

Since  E  (F)=  E  (R)  and  €  =  6  for  a  closed  shell  one 
g  g  v  v 

obtains  E  (R)  -  ©  =  E.(F).  But  E  (R)  -  ©  =  E,(R)  and 
g  v  r  g  v  u 


■ 
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hence  E.(R)  =  E.(F).  Thus  E  (F)-  E,(F)  =6  and  E  (R)-E.(R) 
x<  x-  g  o  v  g  t: 

=  ©v  are  identical  statements  and  the  latter  may  be  taken 
as  the  statement  of  Koopman's  Theorem  for  a  closed  shell  in 
the  formalism  of  the  R  operator. 

The  expression  of  E  -  E^ ,  in  terms  of  ground  state 
quantities  ,  where  either  or  both  of  the  ground  and  ionized 
states  are  open-shell  configurations  can  not  be  shown  in 
general  without  explicit  reference  to  the  spin  states  in¬ 
volved. 


This  can  be  seen  quite  readily  for,  if  the  energies 

are  written  as  in  equation  (4)  of  reference  (18),  terms  such 

as  ff  fV  f^.f^ ,  will  occur.  Expanding  the 

lg  lm  lmg  it  mt  lmt  ^  ° 

integrals  gives 

ig  mgv  img  im  img  im 


it  mt  lmt  im  lmt  lmt 

In  general,  to  express  this  difference  in  terms  of  ground 

state  parameters  only,  it  would  be  necessary  to  relate  the 

ionized  state  a's,  b's  and  f's  to  those  of  the  grourd  state. 

Such  a  general  relation  can  be  simply  found  for  the  f's 

because  they  are  independent  of  spin  states.  On  the  other 

hand  the  relation  of  the  sets  of  a's  and  b's  for  ground  and 

ionized  states  is  highly  dependent  upon  the  spin  character 

of  both  states.  For  example,  compare  the  parameters 

4  2 

required  for  the  description  of  the  S  and  P  states  arising 
from  the  ionization  of  the  common  P  ground  state  of  the 


45 


carbon  atom.  The  relation  between  the  f's  of  the  ground  and 

ionized  states  is  identical  for  these  cases,  whereas  that 

between  the  sets  of  a's  and  b's  cannot  be  simply  generalized, 

Hence,  it  is  suggested  that  a  demonstration  of  the  equality 

Eg  -  E^  =  ,  as  an  expression  of  Koopraan's  theorem,  be 

undertaken  in  the  context  of  the  particular  case  considered. 

One  such  case  which  is  considered  in  Chapter  V  is 

that  of  Li  -  Li+.  Here 
2  2  2 

E  =  22  ,  f .  H.  +  2 .  ,2.  ,f.f  .1.  .  =  2Hn  +  H0+  I, ,  +  I,  0+  I00 

g  i=l  l  l  i=l  j=l  l  j  ij  1  2  11  12  4  22 

1  11 

E  =  22.  ,f.H.  +  2.  , 2  .  , f . f  .1 .  . =  2H,  +  I, ,  so  that 

u  1  -L  1  1  1-LJ-LJ.0  1  <J  1  -L -L 

E  -  E  =  H0  +  I,  0  where  Has  been  set  to  zero. 

g  t  2  12  22 

I? 

^  22  >  =<y,2./K  +  ip  f2>  =  H2  +  I12  and  one 


obtains  E  -  E,  = 

g  t  (v 


IY-3  Virtual  Orbital  Energies 


The  virtual  orbital  energies  for  the  closed  shell 
formalism  of  Roothaan  could  be  used  in  the  calculation  of 
excitation  energies  assuming  the  other  orbitals  are  the  same 
in  the  ground  and  excited  states.  However,  it  can  be  shown 
that  this  application  is  possible  only  for  particular  cases 
in  the  formalism  of  the  general  SC?  treatment. 

In  equations(lII-2,  III-3)  the  k,l  refer  to  occupied 
orbitals  and  the  summation,  of  course,  extends  over  occupied 
orbitals  only  and  in  fact  the  development  of  section  III-l 
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makes  no  mention  of  unoccupied  orbitals.  However,  in  the 
LCAO  context  occupied  and  unoccupied  orbitals  can  be  ob¬ 
tained  and  whether  occupied  or  not  all  must  be  orthogonal. 
Letcing  be  a  non-occupied  orbital  one  has  that 


n  = 

'm 


Z.r£* 
k  k 


y  n 


4 


o 


m 


4 


v/rr><^^x> 


+  ^Cx/kT';! 


Because  of  the  range  of  the  summation,  k^m,  and  by 

virtue  of  orthogonality  { YY °V =  0.  The  equation  reduces 

k.  in 

Multiplying  on  the  left  by  the  bra  <'^m^c'v/  gives 

Again =  0  ,  and  the  right  hand  side  is  zero  giving 

=  o 

But^^W/^^  =  6^*  so  that  the  for  the  unoccupied 

orbitals  have  a  value  zero  simply  on  grounds  of  orthogonality 
and  hence  cannot  be  considered  as  orbital  energies. 

Virtual  orbital  energies,  although  unobtainable  as  the 
eigenvalues  of  the  R*1^  operator,  may  be  found  for  certain 
cases  by  the  use  of  the  operators  which  occur  in  the 

formulation  of  the  . 

For  occupied  orbitals  one  has,  from  equation  (24)  of 


to  R^^=  S 
m 
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reference  (18)  that 

^  y?  yR4  ^/y^  </  y?  /F^<X//>^*J'0^^>  =6^ 


i.e.  the  for  the  occupied  orbital  Vt  ^  "v  can  be  obtained 

n  x  y  l 

with  equal  validity  as  the  expectation  value  of  H*1  ^  or 
But  from  the  definition  of  the  orbital  energy 


F^ 


[i<*  _ 


=  /fT  one  has  that 


F^ 

nf  -of*/  -§-  //f“> 

fi 


Taking  the  definition  of  the  virtual  orbital  energy  to  mean 
the  expectation  value  for  the  virtual  orbital  of  that  ope¬ 
rator  which  may  be  used  to  evaluate  the  occupied  orbital 


energies,  one  obtains 


ft-* 

i 


_  ^  (10< 
mm 


as  the  virtual  orbital  energy.  The  problem  then  is  to  choose 


the 


F^ 

t ' 

In  closed-shell  configurations  all  the  F^^  are  the 


M  0( 

same.  Sven  in  some  open-shell  configurations  the  F7  may  be 
the  same  for  a  given  symmetry  designation  (e.g.  (1S)^(2S)'*' 
(2P)n  ).  The  proper  operator  is  immediately  obtained  in 
both  cases.  However  in  situations  where  the  F^are  dif¬ 
ferent,  the  selection  of  the  proper  operator  cannot  be 
made  unless  some  further  step  is  taken  such  as  using  an 
f(£*  averaged  within  that  symmetry  species. 


•« 


CHAPTER  V 


ATOMIC 


APPLICATIONS 
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V-l  Introduction 

As  indicated  earlier,  properties  of  open-shell  con¬ 
figurations  have  been  obtained  fro™  calculations  upon  the 
nearest  closed-shell  configuration.  The  basic  approximation 
involved  here  is  that  the  orbitals  in  the  open-  and  closed- 
shell  will  not  in  fact  differ  appreciably  and  may  be 
taken  to  be  the  same  (this  will  be  refered  to  as  "constancy 
of  orbitals").  Since  the  formalism  for  open-shell  configu¬ 
rations  is  now  available  it  is  possible  to  test  the 
adequacy  of  this  approximation. 

In  this  chapter  the  approximation  will  be  tested  in 

*J*  *}■ 

the  atomic  systems;  He  ,  He ;  Li  ,  Li ;  Be  ,  Be ,  by  comparing 
the  orbitals  found  for  the  ground  and  the  ionized  states 
and  by  comi>aring  the  derived  quantity,  the  ionization  energy, 
as  given  by  the  approximation  and  by  the  proper  calculation. 

V-2  The  Wave  Functions 

The  total  wave  function  is  given  by  an  antisymmetrized 

r 

product  of  one  electron  functions  _P=  ATl_^  (/L  where  the 
are  given  as  expansions  in  Slater-type  functions  of 
spherical  symmetry 

^(r)  =  2scis^  where  '/Q  =  Hsrs  1  expC-Z^d. 

N  is  the  normalization  constant, 

(2Z  )2s  +  -1 

Ns  =C“'4lTT2iTT~" 


] 


1/2 
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21 

The  orbital  exponent  Z  was  evaluated  using  Slater's  rules 

3 

as  a  guide.  The  values  selected  for  these  systems  are  shown 


in  the  table  below. 

-  Z  - 

System  s 

1 

2 

3 

He 

2.000 

0.575 

0. 100 

Li 

2.700 

0.650 

0.300 

Be 

3.700 

0.975 

0.300 

V~5  The  Integrals  and  the  SCF  Solution 


The  one-electron  integrals*,  H  ,  required,  can  be  formu- 

j/  s 


12 


lated  by  standard  methods  to  give 


H 


rs 


sZ2  +  rZ2  -  (sZ  -  rZ  )2 

■p  o  x  T*  c  ' 

2(r  +  s  -  1)  (r  +  s) 


S' 

rs 


where  s'  =  4ttN  N  [  ^ .  T-S/l  _  ] 

rs  r  s  (z  +  z  }r  +  s  +  1 

r  s 

The  two-electron  integrals*  were  derived  by  successive 
differentiation  of  an  expression  for  the  simplest  integral 

<i-°-<*i(1) Xi(1)/i//1(2)4(2)>  -<& Ks/i-zK, X*> > 


12 


12 


Barnett  and  Goulson  give 


O^kl  ^k2//12//^3  ^KiK2K3K4^ 


5^3xt5>1/2  Tr”3xr“2 


k^k,  *-(  i 

a  b 


-  2  +  ]/>) 


where  the  K 

s 


are  the  orbital  exponents 


Z  and 
s 


*The  one-electron  and  two-electron  integrals  and  a  program 
to  calculate  them  were  developed  by  Mr,  R.E.D.  McClung 
for  another  problem  and  were  adopted  for  the  purposes  of 
this  section. 


- 
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K3  *  K4 
Kb  ■  JL-2-i 


V  = 


*b 


K  "  +"K. 

a  b 


Algebraic  manipulation  of  these  equalities  allows  the 
integral  to  be  cast  in  the  more  convenient  form  as 

<A1(1)/1(1)/  ^2//1(1)^1(1)>  - 

T  r  Ka  +  3Vb  +  Kb  ,2 

x  L — - x~  J 

KaS  <Ka  +  V? 

3  •  _-i  -Z.  r  .  -Z.r 

If  it  is  noted  that  rl  e  1  =  -  r1  e  1  ,  then  the 

integral^  /  Z/2'Z^2^  ^  can  be  derived  from 

°  N/Vp  '  q  r12  r  /s 
the  above  formula  as 

<Z  (D^CD/A.  (1)Z(2)>  =  7t2  2“  (P+A+r+s-4)N  N  N  N 
P  '  0.  '  p.^'  ^  s  /  p  q  r  s 


x 


-jpn 


-2 


^KP+<1 
a 


=2  { 


r+s-2 


[ 


Ka  +  5*a*b  +  4 


K+S^  L  ?T<VV 


A  program*  to  calculate  these  integrals  was  available 
for  the  IBM  1620.  A  master  program  was  constructed  which 
formulated  the  matrix  representation  of  the  pseudo-eigenvalue 
equations  (III-4)  for  both  the  ground  and  ionized  states 
and  solved  the  matrix  equations  in  a  manner  similar  to  that 
of  section  I -7. 

*See  footnote  previous  page. 
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V-4  Results  and  Discussion 


The  calculated  total  energies,  ionization  energy, 
highest  occupied  orbital  energy  and  experimental  ioniza¬ 
tion  energy  are  shown  in  Table  VII.  The  difference 
between  the  ionization  energies  as  given  by  E^  -  and 
by£v,  the  orbital  energy,  is  tabulated  in  the  sixth 
column  of  the  table.  The  percentage  contribution  of  the 
experimental  ionization  energy  to  the  total  energy  is 
given  in  the  last  column.  The  3CF  coefficients  for  the 
ground  and  ionized  state  are  displayed  in  Table  VIII. 

A  comparison  of  the  orbitals  for  the  ground  and  ionized 
states  in  the  three  systems,  shows  quite  clearly  that  the 
constancy  of  the  orbitals  becomes  progressively  better  from 
He  to  Li  to  Be.  The  degree  of  constancy  is  parallel  to 
the  contribution  of  the  electron  to  be  ionized  to  the  total 
energy,  as  displayed  in  the  last  column  of  the  table.  Thus 
in  helium, where  the  orbitals  differ  noticeably , the  ioniza¬ 
tion  energy  is  35%  of  the  total  whereas  in  beryllium, where 
the  orbitals  are  very  similar,  the  ionization  energy  is  2%. 
Hence  it  appears  that  the  approximation  will  become  more 
accurate  the  larger  the  atomic  system  and  in  fact  even  with 
a  system  as  small  as  beryllium  the  approximation  is  already 
quite  good. 

A  comparison  of  the  derived  quantity,  the  ionization 
energy,  indicates  that  the  change  in  its  value  as  given  by 
the  two  methods  is  some  3%  in  helium,  2%  in  lithium  and  1% 
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in  beryllium.  This  trend  is,  of  course,  consistent  with 
the  observations  above  and  indicates  that  the  orbitals, 
even  with  the  small  beryllium  system  are  sufficiently 
constant  that  a  derived  quantity  such  as  the  ionization 
energy  is  given  to  ^uite  good  accuracy  by  the  closed- 
shell  calculation. 


.rtoitsiuolso  lie.'  -’ 
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TABLE  VII 


Constancy  of  the  derived  ionization  energy. 


System 

Electronic 
Energy  (au. 

Ionization  Energies  ( 
)  E  -E.  e  Expt . 

'  g  t  V 

! au.  ) 

E  -E.-e 

g  t  V 

s 

He  + 

-2.000 

-0.787 

-0.903 

0.037 

33 

He 

-2.787 

-0.824 

Li  + 

-7. 224 

-0.194 

-0.198 

0.004 

3 

Li 

-7.418 

-0.198 

Be  + 

-14.247 

-0.309 

-0.342 

0,003 

2 

Be 

-14.556 

-0,312 
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TABLE  VIII 

Coefficients  for  the  occupied  orbitals. 


System 

C2 

C5 

He  + 

1.00 

0,00 

0.00 

He 

0.95 

0.16 

0.01 

Li  + 

1.00 

-0.02 

-0.01 

Li 

1.01 

-0.04 

-0.01 

-0.12 

1.02 

0.00 

Be  + 

1.00 

-0.02 

0.00 

-0.18 

1.03 

-0.04 

Be 

1.00 

-0.01 

0.00 

-0.20 

1.02 

-0.02 

-  . 

CHAPTER  VI 


MOLECULAR  APPLICATIONS 
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VI-1  Introduction 

In  theory  it  would  be  possible  to  carry  out  a  pro¬ 
cedure  exactly  as  in  the  last  chapter  to  test  the  constancy 
of  the  orbitals  for  large  polyatomic  molecules.  In  prac¬ 
tice  this  is  not  so  for  the  complexity  of  the  problem 
forces  the  introduction  of  further  approximations.  Clearly 
it  is  necessary  to  show  that  the  method  of  testing  is 
either  independent  of  these  additional  approximations  or  at 
least  still  reliable  in  spite  of  them. 

Certain  of  the  approximations  of  the  Pariser-Parr 
method  are  investigated  below  aid  numerical  tests  carried 
out  on  the  systems  of  the  previous  chapter  to  determine  the 
implications  of  the  approximations  upon  the  testing  method  , 
Having  established  that  the  test  method  is  reliable  within 
a  context  similar  to  that  of  the  Pariser-Parr  method  an 
investigation  of  the  constancy  of  the  orbitals  is  carried 
out  for  some  molecular  systems. 

VI-2  The  Approximation  of  Orthogonal  Basis  Functions 

In  the  Pariser-Parr  method  approximations  are  made, 
the  effects  of  which  are  difficult  to  investigate  analyti¬ 
cally.  However,  the  approximation  of  orthogonal  basis 
functions  can  be  discussed  to  some  extent  and  although 
justifications ^for  this  assumption  can  be  given, it  can  be 
readily  shown  that  a  true  Hartree-Fock  SCF  LCAO  MO  cal¬ 
culation  is  not  being  performed. 


- 
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A  set  of  non-orthonormal  basis  functions  may  be 

employed  to  obtain  the  integrals  H  ,  J  9  K  and 

b  rs5  pqrs*  pqrs 

thence  the  representation  of  the  Hartree-Fock  operator. 

R  c  =  6  S  c  VI-5 

or, for  the  closed-shell  in  the  equivalent  form  given  by 
Roothaan^, 

F  c  =  £  3  c  VI-6 

If  the  basis  functions  If  are  now  taken  to  be  ortho- 

'  r 

gonal  and  the  same  basic  integrals  retained  in  the  formation 
of  the  representation  of  F  one  has 

F  o'  =  6  c  VI-7 

Equations  (VI-6)  and  (VI~7)  would  be  equivalent  if 
m+F  m  =  F;  nere  m  is  determined  by 

ra+3  m  =  I  VI-8 

Roothaan  has  shown  that  F  is  invariant  only  under  a 
unitary  transformation.  But  equation  (VI-8)  can  be  re¬ 
written  as  3  =  (m+)  ^  =  (m  m+)  \  which  becomes,  if  m 

is  unitary  S  =  (I)  ^  =  I.  Hence  if  m  is  to  be  unitary  then 
3  must  be  a  unit  matrix.  Thus  the  matrix  _F  is  invariant 
under  the  transformation  only  if  the  basis  functions  upon 
which  it  is  based  were  orthogonal.  Thus  in  equation 
(VI-7)  F  is  not  the  representation  of  the  Hartree-Fock 
operator. 

It  should  be  emphasized  that  although  a  Hartree-Fock 
SGF  energy  is  not  found  a  variation  treatment  can  be 
carried  out  and  a  minimum  value  obtained. 
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VI-3  Implications 

The  implications  of  the  above  result  were  investigated, 
by  arbitrarily  setting  S  =  I  in  the  representation  of  the 
He,  He  +  ,  Li,  Li+,  Be,  Be+  systems  as  given  in  section  V~3. 
The  SCF  values  of  the  total  electronic  "energies"  for  the 
ground  and  ionized  states  as  well  as  values  for  the  ioniza¬ 
tion  energy  for  the  various  methods  of  formation  are 
given  in  Table  X.  A  comparison  of  the  orbitals  as  given  in 
Table  IX  shows  that  the  constancy  of  the  orbitals  for  the 
ground  and  ionized  states  is  approximately  the  same  as  for 
the  correct  calculation  of  Chapter  V.  The  quantities  in 
the  last  column  of  Table  X  exhibit  the  same  trend  as  that 
of  column  six  of  Table  VII,  indicating  that  the  magnitude 
of  the  difference  in  the  ionization  energy  should  be  very 
close  to  that  obtained  with  a  proper  calculation.  It  may 
be  concluded  that  although  the  total  energy  is  no  longer 
reliable,  comparison  between  the  ground  and  ionized  states, 
which  is  the  basis  of  the  test  method,  is  still  reliable. 

VI-4  Application 

Since  the  method  of  testing  appears  to  be  reliable, 
even  though  S  is  arbitrarily  set  to  I,  an  investigation 
was  carried  out  to  test  the  constancy  of  the  orbitals  in 
systems  where  the  full  Pariser-Parr  treatment  is  involved. 
The  pyridine,  thiophene  and  isothiazole  molecules  were 
chosen  as  typical  systems. 
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The  molecular  orbitals  were  expanded  in  terms  of 

symmetry  basis  functions  which  could  be  linear  combinations 

of  atomic  2P  orbitals.  For  the  numbering  of  centers  below 
3 


Pyridine  Thiophene  Isothiazole 

one  obtains  basis  functions  for  the  symmetry  designation  A 
and  B  as : 

For  pyridine: 

2*  1/2(2P1z+  2 P5z);/£  =  2'1/2(2P2z  + 

aJ3z;  Xtm 

2-1/2  (2Plz"  22 5z)5  *1"  2'1/2(2P2z  -  2P4z) 

For  thiophene: 


2_1/2(2P2z  + 

2P2Z  ’  2ky 


2PC  )  ;  Xz  =  2_1/2(2P,  +  2P„  ) 


-1/2 


( 2PZ  -  2P„  ) 
v  3z  4z' 


The  integrals  for  these  basis  functions  were  calculated  in 
terras  of  the  atomic  integrals  tabulated  in  Tables  XI  and  XII. 


■ 
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The  values  given  for  these  integrals  were  calculated  from 

24 

data  given  in  the  literature  and  that  these  may  or  may  not 
be  the  best  values  for  the  integrals  is  not  in  question. 

It  is  sufficient  for  our  purpose  that  representations  of  the 
operator  R  be  obtained  which  are  typical  of  those  obtained 
in  Pariser-Parr  calculations. 

The  representation  of  the  R  operator  for  both  the 
ionized  and  parent  ground  state  was  formulated  and  the  pseudo¬ 
eigenvalue  equations  solved  for  the  SCF  solutions  in  a 
manner  similar  to  that  of  section  I“7. 

As  indicated  by  the  results  of  the  previous  chapter 
the  total  energies  may  not  be  meaningful.  Table  XIII  shows 
the  orbitals  of  the  ground  and  ionized  states  and  Table  XIV 

displays  the  ionization  energies  given  by  E  -  E.  and  by 

8 

the  energy,  of  the  highest  occupied  orbital.  Again  it 

is  clear  that  the  orbitals  change  very  little  and  the 
approximation  of  open-shell  quantities  by  those  from  a 
neighbouring  closed  shell  appears  quite  good  for  large 


molecules 


6? 


Table  IX 


Coefficients  for  the  orbitals  in  the 
arbitrary  3  matrix  calculations. 


System 

C2 

C3 

IIe  + 

0. 96 

0.27 

0.00 

He 

0.93 

0.37 

0.00 

Li  + 

0.98 

0.  21 

0.02 

0.97 

0.23 

0.  04 

Li 

-0.24 

0.83 

o.  50 

0.98 

0. 18 

o 

o 

o 

o 

Be 

-0.18 

0.93 

0.24 

0,98 

0.  20 

0.00 

Be 

-0.19 

0.  94- 

0.28 
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Table  X 


Arbitrary  S  matrix  calculations. 


Total  Electronic  Energy  (au.)  Ionization  energy  (  au.  ^ 


ystems 

Calculated 

Experimental 

E  -E, 

s  t 

<c 

V 

He  + 

-2.123 

-2.000 

He 

-2.980 

-2.903 

-.903 

-0.877 

Li  + 

-7.437 

-7.280 

Li 

-7.622 

-7.478 

-0.183 

-0.198 

Be  + 

-14.723 

-14.327 

Be 

-13.000 

-14.669 

-0.275 

-0.271 

E  -S. -€ 
g  t  v 


0.026 


-0.103 


0.004 
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Table  XI 


Two-electron  repulsion  integrals*. 


Pyridine 


10.  5 

7.30 

5.46 

4.90 

5.46 

7.68 

10.  5 

7.30 

5.46 

4.90 

5.60 

10.5 

7.30 

5.46 

4.97 

10.5 

7.30 

5. 60 

10.  5 

7.68 

12.30 

Thiophene 

11.90 

6.58 

4.90 

4.90 

6.  58 

11.08 

7.51 

5.63 

5.39 

11.08 

7.40 

5.63 

11.08 

7.51 

11.08 

Isothiazole 

11.90 

7.09 

4.90 

4.90 

6.  58 

12.74- 

7.81 

5.72 

5. 46 

11.08 

7.40 

5.65 

11.08 

7.51 

11.08 


*The  integrals  are  displayed  in  a  compressed  matrix  array 

I  .  -  fa?.  (1)  2P.(1)  -i-  2P,(2)  2P  .  (2)  dt 

1J  1  1  ri2  J  J 

where  the  subscripts  refer  to  atomic  centers.  The 
integral  values  are  quoted  in  ev.  for  ready  comparison 
to  the  literature. 


. 

1 
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Table  XII 


Core  integrals*. 


Pyridine 

42,02  2.40  --  --  —  2.58 

41.78  2.40 

41.71  2.40 

41.78  2.40 

42.02  2,58 

46.04 


Thiophene 

45c  85  3.00  —  --  3.00 

43.24  2.66 

41,88  2.12 

41.88  2,66 

43.24 


Isothiazole 

46.34  3.00  —  --  5.  CO 

47,50  2.66 

42.19  2.12 

41.98  2.66 

43.31 


*The  core  integrals  are  displayed  as  a  compressed  matrix 

Cid  “  -  W1)  ncore^  2V1) 

where  the  i,^  refer  to  atomic  centers.  For  units  see 
the  footnote  to  the  previous  table. 
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Table  XIII 

Coefficients  for  the  orbitals  of  the  Pariser-Parr  calculations. 


System 

Symmetry 

°1 

°2 

°3 

C4 

C5 

Pyridine 

A 

0.57 

0.52 

0.32 

0.53 

A 

-0.27 

0.49 

0.  58 

-0.59 

B 

0.69 

0.73 

Pyridine 

Ion 

A 

0.57 

0.49 

0.35 

0.57 

A 

-0.25 

0.51 

0.  60 

-0.56 

B 

0.69 

0.72 

Thiopb  ene 

A 

0, 65 

0.61 

0.44 

A 

-0.61 

0.08 

0.79 

B 

0.85 

0.52 

Thiophe  le 
Ion 

A 

0.65 

0.62 

0.44 

A 

-0.58 

0 

© 

0 

0.81 

B 

0.87 

0.50 

Isothiazole 

A 

0.62 

0. 51 

0.32 

0.30 

0.41 

A 

-0.11 

0.63 

0.31 

-0.36 

-0.61 

A 

-0.61 

6.02 

0.53 

0.  58 

0.06 

Isothiazole 

Ion 

A 

0.61 

0. 52 

0.31 

0.30 

0.41 

A 

0.11 

-0.62 

-0.32 

0.37 

0.61 

A 

-0,61 

6,02 

0.  52 

0.59 

0.06 

> 

■?  -?  •  - 
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Table  XIV 

Ionization  energies  (au.)  for  the  Pariser-Parr  approximation 


System 

E  -E, 

s  t 

e  -E.-e. 

s  t 

Pyridine 

-G. 424 

-0.424 

0.000 

Thiopl'e  ne 

-0.394 

-0.397 

0.003 

Isothiazole 

-0.401 

-0.401 

0.000 

- 

. 
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THE  VARIATIONAL  APPROACH  TO  THE  DESCRIPTION 


OF  NMR  AND  EPR  PHENOMENA 


CHAPTER  VII 


IMPLICATIONS  OF  THE  PERTURBATION  APPROACH 
ON  THE  DESCRIPTION  OF  NMR  AND  EPR  PHENOMENA 
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VTI-1  Introduction 

The  structure  of  the  spectra  observed  by  nuclear 
magnetic  resonance  (NMR)  and  electron  paramagnetic 
resonance  (EPR)  may  be  ascribed  to  two  sources.  The 
gross  structure  respresents  the  transitions  of  the  system 
between  two  energy  levels  defined  principally  by  the  inter¬ 
action  of  the  nuclei  and  electrons  with  an  external  field 
which  has  been  modified  by  the  induced  magnetic  field  of 
the  orbital  motion  of  the  electrons.  The  hyperfine  split¬ 
tings  of  this  structure  are  considered  as  arising  from 
the  magnetic  interactions  between  the  nuclei  and  electrons 
of  the  system. 

The  interpretation  of  the  hyperfine  structure,  par- 

22 

ticularly  as  given  by  Ramsey  and  others  for  NMR ,  has 

been  hitherto  quite  successful.  Recently,  however,  experi- 

23 

mental  observations  have  been  made  which  are  in  no  way 

consistent  with  previous  treatments  in  that  they  imply 

negative  coupling  constants.  As  a  result  there  has  been 

considerable  discussion  on  the  shortcomings  of  the  previous 

26 

treatments.  For  example,  McLaughlin  has  pointed  out 
the  implications  of  the  closure  rule  as  used  by  McConnell 

2  o 

Alexander  has  focused  attention  on  the  lack  of  spin- 
correlation  in  Ramsey's  approach  and  Acrivos^  has 
investigated  the  effect  of  the  use  of  molecular  orbitals 
in  place  of  valence  bond  structures  in  obtaining  the  spin- 
density.  However,  all  such  discussions  have  been  within 
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the  context  of  the  perturbation  technique  as  first  applied 
by  Ramsey.  Even  those  attempts,  for  example  the  variational 
approach  of  O'Reilly^,  which  have  departed  from  the  per¬ 
turbation  method  have  been  oriented  towards  providing 
a  description  of  the  phenomena  in  terms  of  concepts  intro¬ 
duced  by  the  perturbation  method. 

This  study  has  investigated  the  implications  of  the 
perturbation  approach  itself.  It  is  found  that  the  nature 
of  the  phenomena  is  considerably  obscured  and  that  a  funda¬ 
mentally  different  approach  is  possible  in  which  electron 
spin  and  nuclear  spin  correlation  is  introduced  auto¬ 
matically  and  the  difficulties  of  the  perturbation  treatment 
pointed  out  in  the  previous  paragraph  avoided.  In  an  il¬ 
lustrative  system,  the  equations  obtained  allow  a  quite 
simple  analysis  of  the  factors  upon  which  the  spin  charac¬ 
ter  of  the  electrons  and  nuclei  are  dependent.  In  a  final 
section  the  relation  of  magnetic  resonance  phenomena  to 
chemical  concepts  ( e . g. , electronegativity)  is  demonstrated 
quite  simply. 


VII-2  The  Hamiltonian 


The  total  hamiltonian,  ,  for  a  molecular  system  in 


29 


the  presence  of  an  external  magnetic  field  ,  can  be  given 


as 


where : 


Ht  =  H1  +  H2  +  H5  +  Hc 


Hl=Sk  TVk  +  IT  Zn  *W  HT  )2  +  V+HLL+HLS+H 


LL  “SS 


kN 
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H  2f-&EkZ:N'<5N  C3(Sk-rkN)(IIJ,rkK)  rjj  (Sk’IN)  TjJ  ] 

H?-  ^  *  VA^hXVIh) 

HC=  W?N  +  Sk*k 


can  be  expanded  as: 

-  H°  +  HLL+  HLS+  HSS+  1313  +  Hlb 

where , 


w- 


1 

rkN 


)  +  V 


eV 

2mc2 


r 


1  1  ^AirX  *  ) 


N" 


"N 


kN 


k'N' 

3 — 


) 


k'N' 


Hlb"  mci  ZWEA  rkN  ^  ZNX  rkN^  °^ko 

V  is  the  interelectron  coulomb  potential,  HTT< ,  H-^g  and 

Hgg  are  the  electron  orbital-orbital,  spin-orbital  and 

spin-spin  interactions  respectively.  H1  and  H1  may  be 

regarded  as  nuclear  spin  interactions  with  the  electron 

orbital  motion  in  that  the  nuclear  magnetic  moment  induces 

electron  currents  in  the  system  which,  in  turn,  set  up  a 

2 

secondary  magnetic  field  experienced  by  other  nuclei.  H 

represents  the  dipole-dipole  interactions  between  the 

nuclear  magnetic  moment  and  electronic  magnetic  moments, 
x 

is  the  Fermi  contact  potential  for  the  interaction 
between  the  nuclear  magnetic  moment  and  the  electron 
spin  density  at  the  nucleus.  The  contribution  of  the  ex¬ 
ternal  magnetic  field  through  the  term  °IN  in  Hc ,  is 
the  potential  energy  of  the  magnetic  moment  of  the  nucleus 


in  the  magnetic  field  about  the  nucleus,  ‘sk  is  the 

analogous  potential  for  the  electron  and  the  external  field,, 

The  terras  in  the  hamiltonian  may  be  grouped  in  two 

classes:  those  which  treat  of  interactions  which  give  rise 

la  lb  2  3 

to  magnetic  reasonance  phenomena  H  ,  H  ,  H  ,  £r  and  those 
which  do  not.  Of  the  latter  is  by  far  the  most  impor¬ 
tant  term.  For  the  former  group  it  has  been  shown  that 
E  is  the  dominant  term  and  is  principally  responsible  for 
the  hyperfine  splittings  ofmagnetic  resonance  phenomena,, 

One  may  then  approximate  as  H^  —  +  H-^  +  Hc  where 

3 

H-pc  represents  the  Fermi  contact  potential,  H  . 

The  relative  magnitude  of  H-pC  allows  the  adoption  of 

the  perturbation  method  in  calculating  its  contribution  to 

22 

the  energy.  Ramsey  first  utilized  this  technique  in 
developing  a  scheme  for  describing  the  energy  levels  ob¬ 
served  in  NMR  spectroscopy  and  latterly  the  perturbation 
technique  has  found  ready  expression  in  the  description  of 
EPR  phenomena.  Subsequently  a  considerable  body  of  material 
both  theoretical  and  interpretive,  has  arisen,  all  within 
the  context  of  the  perturbation  method,  which  purports  to 
describe  NMR  and  EPR  phenomena.  Although  such  efforts  have 
been  fruitful  and  have  been  generally  accepted  they  can  be 
misleading.  In  particular,  the  interpretation  of  the  experi 
mental  data  of  NMR  in  terms  of  the  coupling  constants  of 
Ramsey  has  been  based  upon  concepts  which  can  be  shown  to 
be  the  consequence  of  the  perturbation  technique  and  not 
the  result  of  the  interactions  described  by  the  hamiltonian. 


- 
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VII~5  Implications  of  the  Perturbation  Technique 

Because  the  perturbation  technique  must  be  carried  to 
second  order  to  provide  the  description  of  systems  of 
interest  in  NMR  the  fundamental  nature  of  the  phenomena 
has  been  considerably  obscured,  more  so  than  in  EPH„ 

Hence  a  thorough  discussion  of  the  implications  of  the  use 
of  the  perturbation  technique  in  the  description  of  NMR 
phenonena  is  in  order. 

NMR  spectroscopy  involves  magnetic  resonance  transi¬ 
tions  of  molecules  in  a  ground  electronic  state  which  is 
most  often  a  singlet  state.  If  is  considered  as  a 

perturbation  on  a  singlet  ground  state  wave  function,  /0^> 
then  the  first  and  second  order  perturbation  energies  are 


^0/HFC/i><i/HFC/0> 


Se(1)=<c/hfc/o>  =  o;  £e(2)  = 


where  /0)>,  /l),  .  .  .  . /i/>,  etc . ,  are  electronic  wave  functions 
for  the  system  (i.e.,/i)>  only  contains  electronic  space  and 
spin  variables). 

Substitution  of  the  explicit  form  for  gives 


The  integration  implied  by  the  bra-ket  notation  can  only 
be  over  variables  of  the  bra  and  ket  and  thus  there  will 
be  integration  over  the  electronic  coordinates  only.  Hence 
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the  nuclear  spin  vectors  can  be  removed  from  the  integration 

summation,,  For  rapidly  tumbling  collections  of  molecular 

30 

systems,  averaging  over  all  orientations  gives 

i  \  i  —  -  i  i^-n-  n  ^ 

E 


‘(2)  _y  y  1  j  •  j  y  y  y  ( ion  h  \  2 

LrrN'  3  1N ' ^kll 1/  0  ^  3  ; 

^O/^r^S^iXi/  frkiN.  Si/:L> 


Ei  -  Eo 


¥n,]:n‘  IN,JNN' 


with  J-^,  the  coupling  constant.  The  resulting  form  is 
Ijj*  Ijj,  but  it  is  clear  that  this  form  is  a  result  of  the 
use  of  wave  functions  containing  only  electron  variables 
coupled  with  the  use  of  the  second  order  perturbation  term0 
In  the  Ramsey  approach  the  operator  » ° 

now  takes  the  place  of  H-p^  and  is  used  to  calculate  the 
energy  contribution  due  to  the  Fermi  contact  potential  with 
a  wave  function  which  contains  only  nuclear  variables, 
specif ically, nuclear  spin  variables.  In  the  presence  of 
an  external  magnetic  field  the  perturbation  contributions 
to  the  energy  are  now  given  by  the  two  operators 


tR 


N^N,-L]T  J-N,UNN'  and  H°  "  %  %*NHNIzN 


^FC  ~  I  -^TJ  •  ^ 


where  the  magnetic  field  H  is  in  the  negative  z  direction 
The  wave  functions  upon  which  they  act  are,  for  a  system 
containing  two  nuclei  of  spin  1/2, 

/?!>■  °if°<  +  c12<<S  +  013^  +  °14^ 


A(2> 


=  C, 


oCoC 
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+  c 


22 


?  +  c23^  +  c?zt^etc 


24 
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3  c 

In  order  that  the  representation  of  H-^  and  H 

be  diagonal  the  c.  .  must  be  determined  by  solving  the 

<] 


secular  determinant  |H  -  E  <£l  =  0 

ran  inn 1 

where  Hmn=<VHFC+  HCA^> 

The  results  are  /\^~  (X'cx 

A>>-^ 

A3>=^  (v+po 
A.)  -  tf 


These  /\^  ,  each  of  which  has  a  well  defined  value 

of  >  are  the  simultaneous  nuclear  spin  eigenfunctions 

c  3 

of  H  and  The  observed  transitions  of  NM3  are  then 

I?  vy 

taken  to  be  between  the  /?j  .y  .  Inspection  of  the 

shows  that  the  transition  will  be  accompanied  by  a 

"flipping"  of  a  nucleus  from  *<-+(3  or(?  <k  (i.  e. ,  //^^>  ) 

It  is  clear  that  the  foregoing  description  is  a  result 

of  the  fact  that  the  form  of  allows  simultaneous 

c  3 

eigenfunctions  of  H  and  to  be  determined  which  have 

well  defined  values  of  Z-^I^  and  that  this,  in  turn,  leads 
to  a  description  of  the  transitions  in  terras  of  nuclear 
"flips".  The  development  has  also  made  it  clear  that 
the  concept  of  "nucleus-nucleus  coupling",  as  indicated 
by  the  form  1^*1^,, must  be  more  a  result  of  the  method 
than  of  the  interaction,  for  if  one  considers  the  hamiltonian 
hJjJ  =  Hp  +  Hp,~,  +  Hc  ,  there  are  no  terras  involving  the  product 
of  two  nuclear  magnetic  moments. 
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It  can  also  be  readily  demonstrated  that  the  eigen- 
functions  of  the  harailtonian  v/ill  not  simultaneously 
have  well  defined  eigenvalues  of  2^1^  •  To  see  this  one 
notes  that  in  order  that  two  operators  have  simultaneous 
eigenfunctions  they  must  commute.  Writing  A  =  Hc  +  H-^ 
and  B  =  ^ ,  1^. , 

AB  -BA  =  2N£N,Ek  [(IN*Sk)  IM,Z  Tjt'Z  *Sk  ^ 

%ZN,SkSV(^INVSkVIN'Z  "  IN'ZINVSkV')5V  =  x’y’z# 

=  ”  IN'ZIKXSkX^ 


+  (^INYSkYIN,Z  “  IN'ZINYSkY^ 

=  %%.sk  “in'z;inx^ 

+  SkY('INYIN'z”II^INY^ 


SNSk  ^  SkXINY  +  SkYINX-1 
^  0  >  AB  -  BA  ^  0 

Hence  A  and  B  do  not  have  simultaneous  eigenstates  so  that 
the  Ramsey  perturbation  treatment  has  lead  to  a  concept 
which  describes  the  true  situation  incorrectly. 

An  analysis  of  the  descriptions  invoked  for  EPR  would 
lead  to  a  similar  but  somewhat  less  damaging  judgement. 
In  EPR  unpaired  electrons  are  involved  and  the  doublet 
electronic  spin  state  can  be  successfully  treated  with  the 
first  order  perturbation  term.  As  a  result  the  principal 
imperfection  has  been  that  the  transitions  have  been 
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described  as  "flipping"  of  electrons  -  a  concept  which 
is  invalid  on  the  same  grounds  as  the  ''flipping”  of 

nuclei. 

These  comments  do  not  invalidate  the  perturbation 
approach  but  do  point  out  that  some  of  the  more  common 
(and  one  might  add,  most  useful)  concepts  of  NMR  and  EPR 
have  perhaps  achieved  a  degree  of  acceptance  beyond  that 
necessary  and  that  alternate  approaches  to  the  problem 
could  be  fruitful. 


.=14 


CHAPTER  VIII 


DEVELOPMENT  OF  THE  VARIATIONAL  DESCRIPTION 
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VTII-1  The  Wave  Function 

The  investigation  is  carried  out  within  the  molecu¬ 
lar  orbital  context  where  each  electron  is  considered  as 
being  described  by  an  orbital  function  which  may  extend 
over  the  entire  system  and  a  spin  function  is  associated 
with  each  electron.  The  description  of  each  nucleus  is 
given  by  its  own  particular  spin  function.  The  entire 
system  is  then  described  by  an  anti symmetrized  product  of 
the  electron  orbital,  electron  spin  and  nuclear  spin 
functions. 

The  basic  step  in  the  procedure  proposed  is  the  re¬ 
cognition  that  since  an  electron  or  nucleus  need  not  be 
in  an  eigenstate  of  ,  or  ,then  the  individual  elec¬ 
tron  or  nucleus  may  be  in  some  intermediate  state  which 
is  a  superposition  of  the  possible  eigenstates  of  these 
operators.  This  is  analogous  to  Dirac's  discussion  of  the 
states  of  polarization  of  a  photon  passing  through  a  tour¬ 
maline  crystal.^  Let  /*©*)>  be  the  state  function  dependent 
on  all  the  variables  of  the  system  , where  /  *&}  =  /AjZ^  , 
with  A  an  appropriate  antisymmetrizer .  / /$}  /^)>  , 

where  /  is  the  spin-space  function  of  all  the  electrons 
and  /'Yfy  is  the  spin  function  of  all  the  nuclei.  Let  the 
function  /$ y  be  a  product  of  one-electron  functions 

/0i>-  /h>  /)i>  • 

The  orbital  and  spin  parts  of/0 can  be  given  as  the 
expansion 
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//i>  "  Sil/al>  +  S12/°2>  > 

/6_^  is  a  suitable  basis  set,  /cr^>  and.  /a2/>  are  usual 
spinors,  and  g  ,  for  the  electron.  Thus 


/^^>is  to  be  a  product  of  functions  //£  \  where 


/??p> 


=  E  .h 


d  Pd'  vd 


/Pd) 


where  the  /  ^  . ^>  are  the  spinors  for  the  nuclei  so  that 

d 

/??>=  V%> 


and  hence 


/F>-\/  ^i)VV 


The  nuclei  and  electrons  are  non-equivalent  particles 
and  so  the  antisymmetric  requirements  operate  independently 
on  the  electrons  and  the  nuclei.  Thus  there  will  be  an 
antisymmetrizer  for  the  electrons  completely  independent 
of  that  for  the  nuclei.  The  usual  approximation  that  the 
spatial  coordinates  of  the  nuclei  remain  unchanged  will  be 
invoked^  .  This  in  turn  means  that  whether  the  spatial 
function  of  the  nuclei  is  symmetric  or  antisymmetric  will 
not  matter  to  the  energy.  This  implies  that  no  symmetry 
requirements  can  be  placed  on  the  nuclear  spin  function,  so 
that  either  a  symmetric  or  an  antisymmetric  nuclear  spin 
function  can  be  employed.  Hence  in  this  derivation  A 
will  only  refer  to  the  electronic  antisymmetrizer. 
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7111-2  The  Variational  Equations 


There  are  varying  degrees  to  which  a  variational 
calculation  may  be  carried.  The  most  general  case  will 
be  eriven  by  varying  the  c,  g  and  h.  For  the  total  energy 
given  oy 

w  _oe>/HT/e> 
a  "  ce/e<>  ” 

the  variation  in  the  energy  is  given  by 

£e  =  0  =  <£&/ HT/AiP>  -  E  (ftf/ A^>  VIII-1 

where 


<«/■ 

But  the  variation  in  /0)>is  independent  of  the  varia¬ 
tion  in  /^y  so  that  equation  (VIII-1)  can  be  written  as  two 
equations 

r7  V  s  /-f^/  A^>3  *  o  vin-2 


S^BCfjP  ,/  V  A4^>  -  E  <^-/  a #>>]  =  0  VIII-3 

But  since  0  =  / 0>Y )> 

then 

<tf‘/  ■  <-|f /  +  <^/0 

and  equation  (VIII-2)  can  be  separated  into  two  equations, 
since  the  variations  <f  p  and  S' jj  are  independent.  One  has, 
finally,  .  ^  .  ,T 

A <py\  =0  VIII-4 

£$  a4F>3  =  o  tiii-5 


83 


For  the  ^  ^  of  different  symmetry  equation  VIII-4 
can  be  written  as  n  independent  equations 


/Ht/AJ?>-  E<-^  /  Af>]  =  0  ,  i  -  l,n  VIII-6 


When  certain  of  the  ^  are  of  the  same  symmetry  then  not 
all  variations'^  are  independent  and  a  procedure  similar 
to  that  of  reference  (17)  or  (18)  would  be  required  for  a 
variation  of  the  orbitals. 

Since  no  constraints  are  placed  upon  the  /£$>  or  the 
equations  (VIII-3)  and  (VIII-5)  become 


^k[<L7fr  /VA^-  e<T^77  /Af>]  -  0  ;  i  -  l,n 


VIII-7 


a#->-  s (-4~—/aF>  ]  =  0  ;  P  -  l,q.  VII 1-8 

1  ‘>\-p 

The  method  of  section  1-4  allows  equations  (VIII-6,  “7? 
-8)  to  be  written  in  matrix  form 

0i  ^  0±  _ 

H  1  c  =  E1  S  1  c  i  =  1 ,  n 

H  gJ=  EJ  S  gJ  j  =  l,n 

H  hp  =  Ep  Sp  hp  p  =  l,q 


where  matrix  elements  similar  to  those  of  that  section 
maybe  defined.  This  set  of  matrix  equations  maybe  taken 
as  the  statement  of  the  variational  problem  for  the  cal¬ 
culation  of  all  the  orbitals,  electron  spin  functions,  and 
nuclear  spin  functions.  No  intermediate  operator  (such  as 
Hpg)  is  introduced.  No  restrictions  upon  the  spin  states 
or  spin  properties  of  either  the  electrons  or  nuclei  are 


. 
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implied.  The  orbital  and  spin  functions  are  to  be  deter¬ 
mined  simultaneously  to  give  the  minimum  energy  -  the  best 
overall  description. 

VIII-3  Approximations 

The  development  to  this  point  has  been  towards  a 
quj  te  general  variational  treatment  for  obtaining  the 
lowest  energy  of  H^,  .  However  ,the  observed  energy  differ¬ 
ences  for  the  hyperfine  splittings  are  of  the  order  of 
10  ^au,  Hence  ,in  order  to  obtain  results  of  significance  , 
it  would  be  necessary  to  carry  at  least  16  figures  in  all 
calculations.  Such  a  procedure  is  by  no  means  impossible 
with  the  computational  devices  available  but  a  less  formi¬ 
dable  calculation  will  serve  to  illustrate  the  features 
of  the  variational  approach. 

If  the  first  term,  H^  ,  in  H^  is  dropped,  the  hamil- 
tonian  becomes  simply  a  sum  of  the  Fermi  contact  potential 
and  the  contribution  due  to  the  external  field.  For  an 
investigation  of  the  hyperfine  structure,  only  the  field - 
independent  Fermi  contact  potential  need  be  considered. 

The  implications  of  such  a  rather  drastic  truncation  of 
the  hamiltonian  deserve  comment. 

The  central  point  of  the  variational  method  is  the 
calculation  of  a  minimum  expection  value  for  the  operator 
in  question.  It  is  quite  clear  that  the  minimum  of  H^ 
would  be  determined  chiefly  by  H^  and  so  the  variational 
determination  of  the  orbitals  would  be  principally  the 
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result  of  the  effect  of  H™.  Evidently  a  determination 

Hi 

of  the  orbitals  solely  on  the  basis  of  the  operators 
fip£  and  H  would  be  a  rather  useless  approximation.  For 
this  reason,  when  is  omitted,  as  suggested  in  the  above, 

variation  of  the  is  omitted  and  the  orbitals  used  may 

be  determined  by  alone  by  another  method  (i.e.,  the 
Hartree-Fock  LCAO  MO's).  The  form  of  the  energy  expres¬ 
sion  used  there  will  chiefly  determine  the  spin  functions 
of  the  electrons.  They  can  be  only  slightly  modified  by 
the  Hp,Q.  Hence  if  the  approximation  is  to  be  valid  the 
spin  functions  determined  solely  by  must  be  similar 
to  those  which  would  be  determined  by  +  H^.  This 

can  only  be  done  for  systems  in  which  would  determine 

non-singlet  spin  functions  at  the  minimum,  for  the  minimum 
for  involves  a  non-singlet  spin  function.  Since 

only  influences  the  nuclear  spin  indirectly  through  the 
electron  spin  (the  nuclear  spin  functions  are  determined 
chiefly  by  H-^),  it  appears  that  the  truncation  of  the 
hamiltonian  to  the  single  term  is  reasonable  for  the 

variational  calculation  of  the  spin  functions  of  the 
electrons  and  nuclei  in  systems  which  are  expected  to  have 
non-singlet  ground  states. 


* 


CHAPTER  IX 


ILLUSTRATION  OF  THE  METHOD 
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IX-1  The  Matrix  Representations 

As  an  illustration,  the  system  of  two-electrons  and 
two  protons  was  selected. 

The  orbital  functions  are  then  and  with  the 
electron  spin  functions  ^  and  ^  and  the  nuclear  spin 
functions  and?^,  where: 

^1  =  Si  Ci^  ^2  =  2i  ci' 

^1  =  ^ 2  =  S2  f? 

^  1  =  h^  +  h^  \  2  =  h-[  +  h^  (9 

The  truncated  form  of  the  hamiltonian  discussed  in 
the  previous  chapter  was  used.  The  two  orbitals  ^  and 
^2  are  considered  to  be  fixed  and  a  variational  treatment 
carried  out  only  on  the  electron  and  nuclear  spin  functions. 

The  matrix  equations  to  be  solved  are: 

li_^  Si  ^i_. 

S  6m  -  I  s  gm 


The  total  energy  matrix  E  is  diagonal  and  the  matrices 
H  and  S  have  the  matrix  elements  given  as: 

Hmm  '*^1^2  ara^ 2  Xl ^  2//HFC/  ^  ^1  ^am '  ^  2  "  4  ^2am  ’  ^ vlllr'2> 
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Hnn  •  =<^L ^2  ^lar7ll^2/HPC// ^  ^2^1  an 8  “  ^2 ^Xan '^1^  1^2> 

H^  ^^^l)2(’rVHFc/^1^2^2  “^2^2^1^fr ,7?2> 

Hss'  ^2^1?2\^s/HPc/^  *1^1^  2  “^^J^l^l  ?s  ^ 

Smm' =^4am^2W2/^1^2am,)2  "  ^2^1^2am'  ^1^2 ^ 

^in 8  1^2//^^1^2^ lan'  ”  ^2^1  an')  1^1^2> 

Srr '  ^1  ^2 4^2  ~  ^2^1?2^1')  ^r '^2^ 

Sls  '  ^1^2  " 

The  integrations  implied  by  the  above  are  carried 
out  in  Appendix  IV  and  the  matrix  elements  as  developed 
there  are  displayed  below  as: 

S  matrices: 

S?1 

Sii=(^22  ~  S12slsl')  \ll(22 

S12=  “  S12  sls2^  11^22 

S22=  ^22  "  S12s2g2^  *1  lll22 

Sll  =  ^ll  ”  S12slsl^)??  11^22 
S12=  ~  S12sls2  ^11^22 


89 


0 


si2 

S22 


=  ^22^11  S12  12^  11 


where  Fnm  %s  0<VV  ’  S12  =<^/>2>- 


H  matrices: 


H 


$1 


2  r 


Hn  “  +  DbW  +  (lVD^22Wl 

+  °^1W2]  +  U3CD^22W5  +  DB?11W6  ]  > 

-  2S12(sM[DI\22W5  +  Db\iV  + 

+  Db\iW2])1 


H 


H  .  j£ 


12 


}  22^22*1*  DB?llV  -  S12^l\2\ 


+  DB  ?11W2^  } 

H22  ■  W22CD1?22W5  +  DB?llV  +  (Ul^f?22Wl 

+  +  U3(DI^2W5  +  DbThW6^ 

-  2S12(-6X[d!2722W5  ♦ 

+  db\iw2]  >} 
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H11  -  ■7j-{£ll^DA?22W5  +  DB  !ilWoJ  +  ('lVI,A?22Wl 

+  Db7iiV  +  D6[I)1?722W5  +  DB^lV> 

-  2S12(S^[D^22W5  ♦  D^Wg] 

+  SlS2[Dl2;*22Wl  +  DB2tnW2]  >  } 

Kli  =  ♦  db'?hw2]  -  Sl2tDi2?22ffl 


H22  +  DI?llV  +(U4[Dl722Wl 

+  DB^11W2]  +  U6CD1^2W5  +  DB^11W6]) 

-  2s12(-s”s”  [^2?22w5  +  bI\^61 

+  s“slCI>i2^22W1  +  Db\iW2]  >  } 


H 


?U 


11 


-  4h 


22tDA^22u6  T  vkhlu$ 


1  * 


,U.  +  D 


2^ 


U,  - 


2Dl2si2VlV 


+  W6CDB?22U6+DI^11U3-2DB2r12V1V3]  } 

H12  ■  t%2  [  D^22U4  +  D^11U1  -  2DI2s12V1V2  ] 

*>7  2 

H22  “  \  {’ VD^22V  DI^L1U5-  2Di2£312VlV 

+  'VfeV  DB  ^11U1  -2DB2s12V1V2] 
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+  W6CDbI22U6  +  DB^11°3  "  2D!2si2VlV  ) 

^2  , 

Hll=  T^ll^B^  U6  +  Db4iU3  “  23B  S12V1V3'' 

+  VD1*22U4  +  DA;ilUl  -  2Di2si2VlV2] 

+  VD1?22U6  +  Dl4lU3  -  2D12SX2V1T3]  1- 
Hll=  vKl  DB^22U4  +  DB?11U1  -  2DB2s12V1V2] 

*&’  I  *  -llCD^22U6  +  -  2DB2s12V1V3] 

+  +  Df\iui  -  2Di2si2viV 

+  VD^22U6  +  D^11U3  •  2D12S12V1V3]  > 

Where  the  constants  of  the  problem  are  defined  as: 

dJ  .<J$/£‘(rj)/^>  ;  Dj  -  <-^/r(rB)/^> 

DA  ^ ^  5  DB  -</y^  (rB)/^2> 

D^2=  <^/ r(rA)/^>  ;  Dj2-  <4/  ArB)/f2> 

The  overlap  integrals  are  given  by: 

£■  m  m  ,  m  m 

ai1  Bi  Bp  +  s2  s2 

t  _  n  n  ,  n  n 

>22"  S1  S1  s2  s2 

^12-  V1 


+  h| 
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The  symbols  are  explicitly: 

W1  =  2hl  h2 

W5  =  h[  hj  -  h^  h|  ; 

ui  =  2s?  4  > 

TT  n  n  n  n 

u3  =  Si  Si  -  s2  s2  ; 

Tr  m  n  .  m  n 

vi  =  si  si  +  s2  s2  ; 

v  _  m  n  _  m  n 

oq  Sq  §2  &2  o 

IX- 2  Analysis  of  the  Matrix  Equations 

As  indicated  in  the  introduction  there  has  been  con¬ 
siderable  discussion  of  late  on  the  standard  perturbation 
approaches  to  NMR  and  EPR  energy  levels  in  view  of  the 
disagreement  between  experiment  and  theory,  with  much 
criticism  leveled  at  the  lack  of  spin  correlation,  the 
closure  rule,  etc.  .  The  approach  described  earlier  has 
removed  those  features  associated  directly  with  the  per¬ 
turbation  technique  and  an  analysis  of  the  resulting 
equations  allows  some  useful  discussion  of  the  source 
of  spin  configuration  interactions  (or  spin  correlation) 


W2  =  2h®  h| 


W6 

-  >,s 
"  hl 

■ 

-  v^s 
h2 

h2 

U4 

=  2g^ 

i  m 

s2 

U6 

m 

m 

m 

m 

=  gx 

§i  - 

-  s2 

s2 

TT 

m 

n 

.  m 

n 

V2 

=  Bl 

S2  ' 

*■  §2 

§1 

in  these  phenomena. 

Spin  correlation  (or  spin  configuration  interaction) 
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is  usually  taken  to  imply  that  the  spin  functions  should 

he  given  by  linear  combinations  of  eigenfunctions  of  the 

2  2 

total  spin  operator  S  or  I  ,  for  example,  the  spin  func¬ 
tion  for  two-electrons  would  be  of  the  form 

C1  Ox(l)  (X(2)]  +  c2  [<v(l )(?  (2)  -(?  (1)  <x  (2)  3  +  .... 

■  (g^  +  g2(?)(l)  +  (g,«  +  S 4?)(2),  or,  in  other  words, 
that  there  be  a  ''mixing"  of  the  o<  and  ^  spinors.  If  one 
considers  the  off-diagonal  elements  of  the  H  matrices  as 
effecting  a  mixing  of  the  "states"  indicated  by  diagonal 

^1  5-2 

elements,  then  in  the  H  and  H  matrices  it  is  clear  that 
the  existence  of  off-diagonal  elements  depends  on  or 
being  non-zero.  But  and  can  be  non-zero  only  if  the 
spin  states  of  the  nuclei  are  mixtures  of  the  two  extreme 
spin  states. 

Similarily,  in  the  Hv  and  Hv  matrices  the  mixing 
due  to  off-diagonal  elements  depends  on  and  being 
non-zero,  which  can  only  be  true  if  the  spin  states  of  the 
electrons  are  mixtures  of  v  and  , 

^2 

As  to  the  distinctness  of  the  H  and  W  equations 
(i.e.,  whether^  will  be  different  from^^)  comparison  of 
the  respective  matrices  shows  clearly  that  corresponding 
matrix  elements  will  only  be  different  (except  in  sign, 

possibly)  from  one  equation  to  the  other  if  is  dif- 

2  1  2 
ferent  from  D^,  or  Dg  is  different  from  Dg,  etc.  In 

other  words  the  spin  function  can  only  be  different  if 

the  electron  density  on  at  least  one  nucleus  is  different 


. 
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for  the  two  corresponding  orbital  functions . 

rij  -j  n  o 

Inspection  of  the  H  and  H  matrices  allows  a 

similar  comment  to  be  made  regarding  the  distinctness  of 
and  •  These  spin  functions  can  only  be  different 
if  the  electron  density  associated  with  the  two  nuclei 
is  different. 

In  systems  where  the  electron  density  on  the  nuclei 
is  the  same,  the  result  will  simply  be  that  the  matrix 
equations  are  identical.  This  implies  that  there  will  be 
no  hyperfine  structure*,  a  classic  example  of  this  is  H^. 


IX~5  Sample  Calculation 


31 

Considerable  reference"^  has  been  made  to  the  apparent 
relation  between  the  coupling  constants  for  a  pair  of 
nuclei  and  the  difference  in  electronegativity  of  the  two 
nuclei.  The  perturbation  technique  has  not  lent  itself  to 
elucidating  the  relationship.  It  is  possible,  albeit  in 
a  rather  artificial  system,  to  demonstrate  this  relation¬ 
ship  within  the  context  of  the  approach  introduced  in 
section  IX-1. 


The  D ^  ,  D-g^,etc.,  which  are  constants  of  the  problem, 
represent  the  electron  density  at  the  nuclei  and,  as  such, 
are  a  measure  of  the  electronegativity  of  the  atoms. 

One  then  may  take  the  quantity  ^DA1  +  DA2^  "  (DB1  +  DB2) 

=  ATD^-g  as  a  measure  of  the  electronegativity  difference 
between  centers  A  and  B.  The  "coupling"  constant  of  the 
Ramsey  formulation  is  directly  related  to  the  hyperfine 
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splittings  of  NMR.  Hence  the  relation  of  the  ATD^g  to 
the  hyperfine  splitting  is  equivalent  to  the  relation  of 
the  electronegativity  difference  to  the  coupling  constant. 

The  SCF  solution  of  the  four  matrix  equations  gives 
a  common  minimum  total  energy  E  .  Since  the  basis  set 
is  two  dimensional  a  second,  or  excited  state,  energy 
and  function  may  be  obtained  from  each  of  these  equations. 
It  should  be  emphasized  that  each  of  these  excited  energies 
and  functions  represents  an  excited  state  of  the  system  in 
which  the  other  three  spin-functions  retain  their  ground- 
state  values.  One  then  has  the  energy  level  diagram  of 
figure  7*  Several  arbitrary  sets  of  values  for  electron 
densities  on  the  centers  A  and  B  were  inserted  into  the 
equations  given  in  section  IX- 1  and  the  SCF  energies  and 
functions  were  obtained  from  the  four  matrix  equations. 

The  constants  employed  and  the  energies  and  linear 
coefficients  obtained  are  displayed  in  Tables  XV  and  XVI. 
Figure  8  depicts  the  relation  between  these  electron 
density  differences  at  the  nuclei  and  the  calculated  hyper- 
fine  splittings. 

The  linear  relation  between  ATD^gand  the  hyperfine 
splitting  is  certainly  parallel  to  that  observed  between 
electronegativity  difference  and  the  coupling  constant. 

The  ease  with  which  this  relationship  has  been  demonstrated 
may  serve  to  illustrate  that  various  other  questions  con¬ 
cerning  magnetic  resonance  phenomena  may  be  readily 
answered  when  the  approach  developed  here  is  employed. 
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Table  XV 


Data*  for  illustrating  the  relation  between  ATD^B 

and  hyperfine  splitting. 


db 

4 

da 

da 

1  2 
db  +db 

Di+Df 

tda-tdb 

Splitting 

-  TDb 

-  TDa 

A 

=*tdab 

=  A 

0.3 

o.3 

0.5 

0.5 

0.6 

i— 1 

o 

o 

0.4 

0.8 

0.3 

0.3 

0.6 

0.6 

0.6 

1.2 

0.6 

1.2 

0.3 

0.3 

0.7 

0.7 

0.6 

1.4 

0.8 

1.6 

0.1 

0.1 

0.3 

0.3 

CM 

• 

o 

0.6 

0.4 

0.8 

0.1 

0.1 

0.4 

0.4 

CM 

• 

o 

0.8 

0.6 

1.2 

0.1 

0.1 

0.5 

0.5 

0.2 

1.0 

0.8 

1,6 

*  S]_2  =  0  in  all  cases.  Only  the  relative  magnitude  of 
the  numbers  quoted  above  is  significant  since  no  unit 
is  defined. 
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Table  XVI 


Typical  spin  function"  parameters . 


Function  Parameter  Constants 

1  2 


>1 

0.027 

-0.999 

dA  =  0.3 

*i 

0.999 

0.027 

da  ■  °*5 

OJ 

0.027 

-0.999 

Dg  “0.5 

n 

0.999 

0.027 

d|  -  0.5 

"ii 

0.999 

0.027 

S,£=  0.0 

0.027 

-0.999 

-n2 

0.999 

0.027 

0.027 

-0.999 

*  The  unprimed  functions  refer  to  the  ground  state  and 
the  primed  to  the  excited  state. 
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Figure  7 


Energy  level  diagram 


Legend: 


E  =  minimum  or  ground  state  total  energy. 


The  E 


etc. ,  are  the  total  energy  of  the  system 


when  the  function  of  the  matrix  designated  by 
has  its  excited  form. 


} 


99 


Figure  8 


ATD 


AB 


The  relation  between  the  hyperfine  splitting 
and  the  difference  in  the  electron  density  on  two  centers0 
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EVALUATION  OF  THE  GAUSSIAN  INTEGRALS 
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AI-1  The  Density  Kernel 


To  find  the  density  kernel  for  any  electron  it  is  first 
necessary  to  proceed  through  some  extremely  cumbersome  nota¬ 
tion  and  untidy  equations  for  one  center  with  many  electrons. 
Hence,  the  method  will  be  indicated  by  finding  the  density 
kernel  for  electron  one  when  there  are  two  electrons  on 
one  center.  Although  , as  noted,  the  notation  and  form  of 
equations  can  not  readily  be  extended  to  many  electrons,  the 
procedure  and  sequence  of  steps  which  would  need  to  be  car¬ 
ried  out  should  be  evident. 

The  density  kernel  for  electron  one  for  the  system  of 
two  electrons  on  one  center  is  given  as: 

V  <01/R(r',  ~r"1)/0^)  =  JexpC“B“ir,2"mjr|  ”  ns^'  “  *2 1  2  3 

exp[-nK  r"2-  m_.r2  ~  ns , | r"-r2|  ^Idr^ 


= jexp[-r ' 2(mi+  ng)  -  r"2(mi,+  ngl)] 


exp[-r2  (nij  +  ,  +  ng  +  ng ,  )  +  2fx2<nsr'  +  ng  ,r")]  dr^ 

The  integral  part  may  be  written: 

1 2=  (expC-S^g^r^  +  2r^d^_]  dr2  where 


S1  =  mj  +  ns 

d,  =  n  r' 

1  s 


=  +  ns’ 

d. 


■2  "  V" 


2  2 

r  P  d,  df:  df; 

l2=Jexp[-2tgt(r2-  2?2  “  +  ^ ~  p  3  dr. 

u 
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dt  .  d2 


=  exp[Et  ] Jexp[-Etgt(r|  -  2%  ^  +  -§)]  dr"2 


Sh 


^•4-  /"  ^  P  _v  — *  (“*'+- 

I2=  exp[E  ]  J  exp[-Etgt  (r2-  C^)"]  dr2  where  Qt=  — 
t  t 

This  integral  can  be  evaluated  as  follows: 

I, =  (dr  exp[-E  q  (r  -  Q  )2  where  Q  is  either  some  variable 
1  ;  m  ^  s^s  v  m  ^s  & 

not  dependent  on  r^  or  is  a  constant. 

T-i  =  (dr  exp[-  E  q  (r2+  Q2  -  2r*Q  )  ] 

1  J  m  ^  sus  m  *s  m^s' 


=  f' 


drm  exp[-  Vsrm  -  22^^^)  r^Qs  + 


u  Mu. 


u 


“  S*3(SUtQt)2  +  ] 


UuU 


/*  p  ^t^-t  — *■ 

=  (dr  exp[-  E  q  r  -  2E  q  — ^ —  r-Q. 
J  m  ^  sHs  m  s^s  Eq^  m^t 


+  S^s^S. 


Ssqs\  -,2 


q 

u^u 


)  -  2oq  ( 


sHsvE  q  )  +  ^s^-s^s  ^ 


_ r  r.  „  ^_2  0  /^t^t^t  ^  t  -j 

)drm  exp[-  Ssqs(rm  -  2  (-g-r - ).  +  (— — ~)  3 


.  q  '*  m  '  E  q 
u^-u  uHu 


exp[ 


Ssqss(EtqlA) 


s  q  S  a  SsqsQs  ^ 


U"U  UAU 


=  (dr  exp[-s  q  (r  -  )2  -  E  q  Q2  +  (^3) 2  ] 

j  m  ^  s^s  m  E  q  '  s^s^s  v  E  q  '  J 

^u  u^u 

=  exp[-E  a  Q2  +  C-—^)2]  (d?  exp[-E  a  (rm- 


s x s  “s  '2udu  '  “  >  m 


s  s  m  Vu 


By  this  sequence  of  steps  the  integral  has  become  one  of 
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K 


the  form  (dr  exp[-(r  -  K)^  E  q  ]  since  the  portion  within 
m  m  s  s 

the  bracket  no  longer  depends  on  s9  i0e0  K  is  a  'constant 

=  g"  ;  so  that,  [drm  exp[-b(rra~  K)2],  where  b=Egqg  . 

For  this  integral  it  is  noted  that  as  the  vector  r  sweeps 
out  its  volume  element  the  vector  r  -  K  also  sweeps  out  its 
volume  element.  Hence  the  integrals  may  be  equated  as 
MnkexP["brmk]  “  fd^expC-br^] 

which  has  the  value  (— ?-)^2  =  (^—  )^^2  .  This  in  turn 

qs 

means  that  the  total  integral  1^  is  equal  to 

V  +  (2%.  )3/2 


-  q 

u^u 


S^S 


This  can  be  written  in  a  slightly  different  form  if  it  is 
noted  that 

(2t^tV2=  A^tAA^tA  =  AqsAAqtQt)=  2st^sqtQA 

and 

Sq  Q2=  —  E  q  Q2  =  Eq  Q2 

^■s  s  E.  q  s^s  ^s  Eq  .  ^s^s 


u  xu 


Lu 


=2^(SsEt<ltqsQf  +  SsVsq<)=2ZF  SsSt<1s'1t(Qs  +  '5t) 
So  that  the  integral  1^  may  be  written 


Il=(^-)5/2exp[ 


E  q 
u^u 


T“SsVsqt(Qs  +  V  +  SsVsqM  ] 


^  Sq  ; 


V2expC_  K h  ssWt  («2+  ^  -  2%v 


] 
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-  _1Z2 


=  exp[ 

s^s  u^u 


STq,;  ^sMs^t  (Qs~  V  ] 


■  (rV)5/2  expCVt  qs\(!s~ 

udu 


Finally  one  then  has 

Jdrraexp[-Ssis  (rm-  Qg)2]^^  )3/2exp[S2tqsqpQs-Qp2/:Cqu] 


S^S 


Substitution  of  this  form  of  the  integral  into  the  last 


equation  for  gives 
d2 

V  exp[£(^)]  (g-|-)5/2  exp[iSsqsqt(qt-  Qs)2Atlt] 


t&t 


.  n^/2expfpsr '  Km.i+  V4  (nar")  (mi'+  ns-)_1 


(m  .  +  m  . . +  n  +  n 
3  3  s  s'y 


-] 


n  r 
s 


n  ,r 
s 


(m.+  n  )(m..+  n  .  m.+  n  m  .|+  n  , 

expr-  — J _ s  ,j _ §_ _ a _ s _ i _ §_i 

2L  (m  .+  m  . ,  +  n  +  n  , ) 

3  3  s  s' 


And  the  total  density  kernel  for  electron  one  is: 

3/2 


n- 


2 

jj2  exp[-r 1  (mi+  ng)  -  r"  (mi,+  ng , ) 


-  -r."2/ 


(m  .+  m  . .  +  n  +n  ,) 
a  3  s  s' ' 


n  r  *  nr" 

2  p  s  _  s 

(n  r1)  (n  r")  (m.+  n  )(m.,+  n  .)  Lm.+n  m.,+  n  , 

s  y  +  s  y  _  ,i  sx  3  s 1  y  ,j  s  ,1  s 1 


m  .  +  n 
3  s 


m  .  +  n  . 
3  s 1 


m  .  +  m  . ,  +  n  +  n  , 
3  l  s  s 


2 

■] 


=(  _ 2 _  )3/2  expr-r'2(m  +  "j”3  ) 

'  m  .  +  m.,+  n  +  n  ,  ^  i  m.+  n  ; 

- '  is 


3 

--r."  ^ 


"S  s 

m  •  ,n  , 


r,^(m  +  innsi  ).  ,(y  V  . 

'■  1  V  V  ns  ■ )  l-mj +ns 


ns'r 


m  .,+n  , 

3  s' 


It  is  advantageous  to  have  the  density  kernel  in  the  form 
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Ip  -  G’i  exp[-  ar'2  -  br"2  +  c(rl  “  r")^] 

The  required  values  of  a,  b,  c  and  G  may  be  obtained  by 

expanding  the  square  in  this  last  equation  and  the  square 

in  the  one  previous  to  it,  collecting  the  coefficients  of 

r'  ,  r"  and  r^r",  equating  the  appropriate  coefficients 

and  solving  the  simultaneous  equations  to  give 

n  (n  +  n  . ) 

al=  mi+  ns  "Cm'Vrn  VV°Es-“^T 

J  J  s  s 
n  , (n  , +  n  ) 

bl"  mi'+  ns'  (m^  +  m_. ,+  ns+  nG. ) 


C1  (in  .  +  m  . ,  +  n  +  n  T) 

G  O'  s  s1' 

G,  =  (m  .+  m  . ,  +  n  +  n  ,) 

1  j  J  1  S  s  1  ' 

It  should  be  emphasized  that  although  other  rather 
obvious  approaches  to  the  solution  of  the  integral  may  be 
used  in  this  simple  problem  the  particular  approach  in¬ 
dicated  is  the  most  satisfactory  in  the  more  complex  case 
of  many  electrons  and  many  centers. 


AI-2  The  Potential  Energy  Integral  for  Electron  One 

I3  ■  <01/l/r1/02) 

And  from  the  definition  of  the  density  kernel  this  becomes 

Iy  (dr  P(r/r')  P(f/r")[y«]  <'01/R(r '  ,5"  ,l)/02>  . 

For  the  above  form  of  the  density  kernel  this  becomes 

I^=  G1  jdr  P(r/r '  )  P(r/r")[p.]  expOap'2-  bp1"2-  c(rj-  rp2] 

Q. 

=  G1  [dr  i  exp[-r2(a1+  bx)>  2it 


109 


AI-5  The  Kinetic  Energy  Integral 


The  kinetic  energy  integral  for  electron  one  is 

2 

I4  =  ^1/"  “Tsi  ^“^T’^2^2/>  where  2i  is  the  Lap" 

lacian  operator  in  cartesian  coordinates  r  ,r  ,r  ,  and  in 

terms  of  the  density  kernel^  1^  hecomes 
G 

I,.  = 


Sxyz  ,fd?  P(?/r')  P(?/f")  (-|^)2 

X 

'  2  -  b-^r"2  -  c^(r'  -  r")2] 

-  jdrKr/r' )  P(?/r'')(-|p,)(-2br£  -2cr"  +  2cM 


exp[-a^r,t_  -  b^r' 


exp[-ar'2-  br"2-  c  r ' 2+  r"2-  2(r'r"  +  r'r"  +  r'r")  ] 


x  x 


-  G  Sxyz  jdr  P(r/r ' )  P(r/r") 


y  y 


Z  Z' 


exp[-ar '  2-br,,2-c(r '-r")2]  -[  b+c-2  (brM+  cr"  -  cr')21- 

L  X  X  "X  -J 

Operation  of  the  replacement  operators  may  now  be 
carried  out  to  give 

I4=  Gl2i-fd^  exp[-r2(a+b)] -[  b1+c1-2b2(rx)2J- 
=  G1  jdr  (3b1+  3c1~  2b2r2)  exp[-r2(a1+  b-^  ] 


=  C3b1+  3c1+  2b2  -§* 


]  Jdr  exp[-r2(a1+  b1)  ] 


=  Gl[3b1+  3c1+  2bx  ^ 


2  5  ~\  f  n _ x3/2 


]  ( 


ai+  V 


"  Gi  ( 

G 


71 


I'M  m05/2(3V  3ex)  -  5t>2 


TU 


3 


1  M1 
3/2 


(ai+  V 


5/2 


(a-^  bx) 


372  3(axV  Vl  +  alCl} 
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AI-4  Electronic  Repulsion  Integral 

The  repulsion  integral  for  electrons  one  and  two  on 
the  same  center  has  the  form: 


2  2  2 

I^=  <^exp[~miri  -  m^r2  -  nsr12] 


1 


id 


exp  if  -  mj.rf  -  n^ifd,]  > 

-  (dp  exp[-Urf]  jdr^  ——  exp[-?r|  - 


Wr12] 


where  U  =  m.  +  m.  , 
l  i 

V  =  m  .  +  m  . , 

W  =  n  +  n  . 
s  s ' 


The  inner  integral  is  defined  as 


I5i=  fd? 


2  exP  [_Vr2  "  Wr12] 


and  it  is  noted  that  for  this  integration  electron  one  is 
fixed. 

From  the  law  of  cosines  one  has 
2  2  2 

r2=  r12  +  ri  “  ^r12rl  cos  0>  w^ere  0  is  ihe  angle  between 


r2  and  r,£  • 


It  is  also  seen  (as  for  1^)  that  as  r2  sweeps  out  its 
volume  elements  that  r^2  also  moves  through  its  volume 


element . 


r 


ic.  = 
5i 


exp[-Vr2  -  Wr-^2]  dr2  cirl2  exP^-^r2  ”  ^r12^ 


12 


12 


lr— -■  sin  6  d©  d 0  dr^2  exP^“Vr2  -  Wr^2] 


Ill 


C  2  2 

=  2u  J  sin  ©  d©  cirl2  exPE“Vr2  ~  Wr12^ 

2 

On  using  the  cosine  law  form  for  ^  this  becomes 
Irj,=  27U  |exp[-Vr2]  exp[-r^2  (V  +  W)  +  V(2r^2r^cos 


rl2^r12  s^-n®  ^ 
Integrating  over  ©  gives 


I,..- 


exp[-Vr^] 


5i  Vr. 


rv 


rjf  dr12  exp[-(V+W)  r*,  +  2Vr1r2]-jdr12  ^ 


exp[-(V+W)  r^2  ”  2Vrlrl2^ 

On  completing  the  squares  in  each  exponent  in  the  integral 

this  becomes 


In,= 


expOVrp 


Vr. 


Vrl  1 

expCv+ 

Vr2  / 

exP[V+W]  J 

dr12 

VWr2 

/-00 

exPC  y+W 

r._ 

/d 


o O 


Vr. 


exp[-( V+W)y ' 2]dy ' 


-a 


with  y  =  r-^2 


Vr. 


Vr. 


Vr. 


V+W 


;  y' *  (rq2  v+w  )  ;  a  =  v>W" 


ic.= 

5i 


VWr^ 

7iexp[__V+W  ^ 

Vr. 


oo 


exp[-(V+W)y2]  dy  +  J exp[-(V+W)y2dy 


-a 


x'00  ^co 

j exp[-(V+W)y 1 2dy '  -  j  exp[-(V+W)y '  2dy ' 


*00 


oO 


The  integral  j  exp[-( V+W)y^dy  =  | exp[-(V+W)y ' 2dy ' 


o 


o 
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since  the  y  is  really  a  dummy  variable  so  that  the  value  of 
the  integral  depends  only  on  the  limits.  Hence  the  last 
line  can  be  written  as 


I,.  = 


TC 


VWrJ 
expC”v+w  -1 


5i  Vr 


11 


o 

Jexp[-(V+W)y '  3dy  +  J^exp[-(V+W)y?dyL 
.o'7  —  a  J 


y'  =  r^2+  a  and  for  the  allowed  range  of  the  first 
integral  C^y'^a;  therefore  0><^r^p+  a  ^a  and  hence 
-a  a  *  ^=r12”  a  an<^  ^or  ^e  a^owec^  value  of  the 

second  integral  -a^y^O,  therefore  -a^r^a^O  and  hence 

_a^r12^a- 

Thus  in  the  two  integrals  has  the  range  -a^Tr^o'^-3 

(i.e.  ,the  ranges  of  y  and  y'  have  the  same  range  for  r^2). 
The  two  integrals  can  be  combined  in  terms  of  the 


variable  r^2  as 


I  = 

5i 


TZ 


YWr£ 

6XP[  y+yy  ]  ( 

~l  -a 


exp[-(V+W)r212]  dri2 


which,  since  exp[-( V+W) r^2]is  an  even  function  , gives 

P_  VWr2  f:  P 

I5i=  W~~  exP["  T+W“]  j  exp[-(V+W)r^2]  drl2 
^  1  o 

V 


With  the  substitution  of  r 
becomes 


12  "  v  +  w  z>  the  integral 

2 


271  _  r  _  VW  „2t  r  1  „  r  _  ^  _2n 

r-j^CV+W)  exp^-  rl^  j  exp^-  V+W  z  ^  dz 


The  total  integral  on  substitution  for  becomes 


I5=i'd?lexpC-Url]  §7JJW7 


J^1exp[- 


V2  2-, 

V+W  z  ^  dz 


exp[-  r2n 

t?XpL  y+yy  X]_-> 
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2lL_  Cjf  1_  exp[-Ur2  -  r2  ] 

V+W  Jarl  rx  expL  url  V+W  rl  J 


1  expO  ?fw  ]  dz 


v+w 


dr  1 
drl— r 


rl 

2  r  x  2 

exp[-Gr^]  ^  exp  [-Hz  ]  dz  with 


(UV  +  UW  +  VW)(V+W)  =  G  and  H  = 


r 


V+W 


V  WW  57  expC'Grl]  f  expC-Hz2]  dz 


/  V  s2  2tc 

^  v+w; 


V2 

(V+W) 


^dr-^  —  exp[-Gr^]  ^  exp[-Hz2]d: 


«(JL)2  r_22L 

^V+W;  L  H 


jdr^  ~  exp[-Gr2]  j  exp [-Hz2]  ] 


The  integral  j  dr^  ^  exp-[Gr2] 

12  u 

found  by  standard  methods  to  give 

5/2  .  2U)5/2(p+V)-l/2 


exp[-Hz  ]dz  is  readily 


I  =(^-)2 


2(n) 


5  'V+W'  <JH(G+H)1/2  "TUV  +  WV  +  U’,») 

as  the  form  of  the  electron  repulsion  integral 


which  is  taken 


AI~5  Overlap  Integral 


The 


overlap  integral  1^  is  given  for  two  electrons  and 


the  method  is  readily  extended  to  more  than  two  electrons 


I6=  <(01/02]>  =  CT(  J~  dr\  )  exp[-mir2  -  m..r| 


’  nsr12] 


exp [-UK  , r^  -  ra.,r2  -  ns,r^2  ] 
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exp[-£r^a^ 


where 


£ 

id 


ai = 

(mi  +  mi, )  ; 

a2  ' 

(m  .  +  m  .  . ) 

0  0 

(n  +  n  ,  ) 
s  s 

(n  +  n  , ) 
v  s  s '  ' 

c12- 

2 

c2l” 

2 

2 

1  =  C1T 

6  i 

expC-ZPr^  (ai+ 

2S.c.  .)  + 
i  10 

2Zid(rirP  0 

2  r 

=  ('TT  dr.)  exp[-£.  .r.r.d.  .] 
±  J  i  ^  10  i  J  10 

where 


iP 


d..=  a .  +  =E .  .c  .  . 
ii  i  ij  10 


d .  .  =  -2c  .  . 
id  10 


I  =  1J  r IT.  (dx.  )  exp [-2.  .x.d.  .x  .] 
b  xyz  iJ  i'  1  ij  i  10  0 

To  perform  the  integration  of  the  last  line  it  is 

first  necessary  to  show  that 

1/2 

^x.d.  .x  .1  =(Ti.  I  dv.  )  exof -Evfc  .  1  =  — 


(1TifdxP  expC-S  Xid  x JKTTifdyp  expC-Eyfc^  - 


1Ti°i 


This  can  be  seen  as  follows:  the  integral  may  first  be 
written  in  terms  of  the  matrix  D  and  the  vectors  x  (whose 
components  are  x^,etc.)as 

(Tlj"  dx^  exp[-2ijxidi^x^]  =  (Di|'dxi)  exp[-x+Dx]  . 

Patting  x"  =  •  x+=  y+£+and  therefore  x+_T  x  =  y+£*J^  jj>y=y+ By" 

9  j  ^  ? 

with  B  =  £+D 

Clearly  if  B  were  diagonal  then  y+By  would  be  a  sum  of 

squares  and  hence  the  desired  form.  The  problem  then  is  to 

obtain  the  which  will  diagonalize  the  matrix  B.  £  can 

13 

be  obtained  by  a  standard  method  . 
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For  the  case  of  n=2 


D  -I*11  Dl2|  so  that  2  = 

lD21  °2?J  V° 


1  ”  ^i2//,^lli  +  f3n 

lf-  111  and  £  D£=|  11 


therefore  x+Dx  =  y+j£+D  _^y=  y+  ^ 


D 


11 


0 


=  ^iBii^i  where  Bn  =  DH»  B22=  D 


bl 


11 


0 

di  )  y 


D 


11 


) 


0 

Id 


D 


11 


Hence 


('rlij'dxi)  exp[-S.  -  (liifdyi)  expC-S.B^y2] 


n 


1/2 


71 


1/2 


CTi  b  .  .) 

1 1  ii; 


172 


I  iBii  ^Dll^  ^D22  D12')  =  D11D22  3 

Then  the  integral 

Wz  exP["SijxididXd] 


12 


W 2 

=  IdI  . 


since 


n  1/2  ^ 

-  (TTT72)  n 


71 


5/2 


Dl^  (|D|; 


~W 


for  n  =  3 


Evaluating  this  for  two  electrons  one  has 


Ir  =  i-- 


n 


•1 


5/2_ 


TC 


6  Dll322  D12D12 


( m .  +m .  ,’+n  +n  ,  )  (m  .  +m  .  ,  +n  In  ;-(n  +n  7) 

v  i  i  '  s  s'/vj  j  '  s  s ' y  s  s '  ' 


T  -(  11  ^5/2 

±6"VUV  +  (U+V)  W  J 

Which  is  the  final  form  of  the  overlap  integral. 

In  Boys'  paper  he  requires  that  D  be  a  positive  definite 
matrix  in  order  that  a  real  number  be  obtained  by  the  square 


3/ 
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root  operation.  This  is  true  if  the  parameters  forming 
U,V,W  are  all  positive  (see  denominator  in  the  last  line 
above) . 


APPENDIX  II 


SOLUTION  OF  THE  MATRIX  EQUATIONS 
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AII-1  Solution  of  the  Matrix  Equations 

The  solution  by  the  Jacobian  technique,  of  the  eigen¬ 
value  matrix  equation,  Hcf  =  ESc,  where  the  basis  functions 
are  not  mutually  orthonormal,  must  be  approached  by  trans- 
forming  this  equation  to  the  equivalent  form  H'X  =  EX 
where  the  basis  functions  are  orthonormal. 

The  relation  of  the  two  equations  can  be  established 
in  the  following  manner.  Defining  a  matrix  T  by  T+T  =  S 
and  substituting  in  He  =  ESc"  gives  :  He"  =  ET+Tc.  Inserting 
I  =  T  ^T  into  the  equation  leads  to  HT  "'"Tc"  =  ET^Tc".  This 
form  reduces  to  (T  +  )  ^HT  ^Tc^  =  ETc  on  left  multiplication 
by  (T  +  )  1.  Defining  m  =  T  ''"allows  the  equation  to  be 
written  as 

H'X  =  EX  with  H'  =  m+Hm  and  X  =  m  ^c 

In  asking  for  the  solutions  of  H'c“  =  ESc  one  is  asking 
for  the  eigenvectors  and  eigenvalues  of  the  set  of  equations 
He =  E  S  C^ 

He ^  =  E^ScV*  ,etc. 

a  set  of  equations  which  can  be  represented  by  the  one 
matrix  equation  HC  =  ESC  where  0  is  now  a  square  matrix  and 
E  a  diagonal  matrix.  The  Schmidt  Orthonormalization  method 
(SOMS)  given  in  the  next  section  can  then  be  employed  to 
obtain  the  m  matrix  required  to  transform  this  equation 
to  the  form  H'A.  =  EX. 

The  Jacobian  Technique  described  in  section  AII-3  gene¬ 
rates  a  X  such  that  X+H'X  =  H"  is  a  diagonal  matrix  with 
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A  A  =1. 

e^=> 

Writing  out  the  equation  k  H “A  =  EX  X  element  by 
element  gives 

t^H^i=  E1 
"^2-^2=  S2 

so  that  the  vectors  are  eigenvectors  of  the  matrix  _H' 
with  eigenvalues  E^.  But  these  last  equations  become 
A.m  HmA . =  E.  on  substituting  for  H'.  Writing  c.  as 

i - li  °  —  °  l 

mA^  the  equation  becomes 

c+Hc.  =  E. 

l—i  l 

and  the  vectors  cA  are  the  eigenvectors  of  the  matrix  H 
with  eigenvalues  E^„ 

Since  it  is  the  lowest  eigenvalue  that  is  required  a 
search  (EVSS)  of  the  eigenvalues  E^  must  be  carried  out. 

The  vector  c  corresponding  to  the  lowest  E  will  be  given  by 

«»N 

mA^  where  A^  is  the  row  of  the  A  matrix  corresponding  to 
the  row  in  E  in  which  the  lowest  E^  occurs. 

AII-2  Schmidt  Orthonormalization  (SOMS) 

The  orthonormal  set  of  functions  generated  by  the 
Schmidt  Orthonormalization  procedure  from  a  non-orthonormal 
set  p/  are  related  for  a  set  of  n  functions  by 


V?  =  L  .  .m.  Sfi'  i=  1.  .  . 

_L  r  J  -L  -L 


n 


AII-1 


n.  f  ]°  Nj 

Such  a  form  introduces  — — L  parameters  m..  .to  be 


"2"  *  |  10 

determined  subject  to  the  7 orthonormality  conditions 
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<^/yV  fi3. 

Expanding  gives 

a  «J 

=<^k  mik^>  k  /“I 


=  where  sf  -Of//’.} 

Taking <(y?/  =  as  a  matrix  element  one  obtains 

_S  =  m+S'm  and  the  matrix  m  has  the  triangular  form  required 
by  equation  AII-1. 

In  view  of  this,  the  matrix  operation  m+S'm  to  give  a  unit 
matrix  for  the  solution  of  the  matrix  equations  can 
indeed  be  taken  as  simply  a  transformation  to  an  ortho- 
normal  basis  set. 

The  explicit  form  of  the  matrix  element  m. .  is  most 

conveniently  found  by  first  finding  T . . ,  such  that  T  T=S, 

J 

for  simple  cases  and  then  obtaining  by  induction: 


T.  =  (S.  -  2.  . T,  .T,  ) 

m  v  m  k<i  ki  kny 


1/2. 


'T.  . 

li 


T  =  (S  -  2,  T2  )1/2 
nn  nn  k<n  kn' 


Since  _m  =  _T_  ,  one  can  obtain  the  elements  of  m  by  the 

identity  mj?  =  I;  obtaining  finally  rrr^  =  1/T^ 


m  =  -Sn  1  m . ,  T,  /T 
(t<n)  k<i  lk  nn  • 


AII~3  Jacobi  Diagonalization  of  H1 


The  problem  is  essentially  that  of  finding  the  matrix 
\  which  will  reduce  the  off-diagonal  elements  to  zero  or  as 
close  to  zero  as  required  in  A.+H'\  =  H"  . 
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The  largest  off-diagonal  element  is  determined  and  the 
matrix  found  which  will  reduce  the  value  of  this  element 
in  to  zero  with  the  orthogonal  transformation  =H^  . 

Although  this  transformation  certainly  alters  the  values  of 
the  other  elements  in  it  is  convergent  for  the  change  is 
primarily  towards  an  increase  in  weight  along  the  diagonal 

14 

and  in  fact  it  can  he  shown  that  for  any  off-diagonal 
element  S'rs  .  The  process  is  repeated , each  time  the 

largest  element  in  the  transformed  matrix  being  determined 
and  a  new  \m  found  until  finally  the  largest  off-diagonal 
element  is  effectively  zero  and  the  matrix  X  =  ^2.— 2***— n 
has  been  found  which  diagonalizes  H'  as  \+H'\  =  HM ,, 

The  formulation  of  the  matrix  \  is  given  explicitly 
by  G-reenstadtP'^and  a  method  developed  there  for  solving 
the  trigonometric  equations  which  is  stable  against 
singularities  and  rounding  errors  of  machine  calculations. 


APPENDIX  III 


VARIATIONS  OP  THE  BASIS  SET 
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AIII-1  Mavle 

The  logical  question  of  the  direction  of  change  for 
a  better  value  of  a  particular  exponential  parameter  m  is 
extablished  and  answered  by  changing  m  by  -  Am  and  cal¬ 
culating  the  energies.  Knowing  this,  change  is  carried 
out  in  the  x^roper  direction  until  the  energy  begins  to 
rise.  At  this  point  the  region  of  the  minimum  has  been 
established  and  a  value  for  m  on  either  side  of  the  mini¬ 
mum  determined. 

The  value  of  the  increment  is  then  decreased  by  a 
factor  of  two  and  the  proper  direction  of  change  again  es¬ 
tablished.  As  before,  alteration  of  the  parameter  is 
continued  in  the  proper  direction  until  the  minimum  is 
bracketed  -  this  time  more  closely  than  before.  The  pro¬ 
cess  is  continued  until  the  value  of  m  for  the  minimum  has 
been  determined  to  any  desired  accuracy. 

Having  established  a  temporary  best  value  for  this 
parameter  a  search  for  a  best  value  of  the  next  parameter 
is  initiated  and  the  value  found.  'The  procedure  is  con¬ 
tinued  until  the  temporary  best  values  for  all  the 
parameters  have  been  determined.  A  new  cycle  is  then  begun 
with  the  variation  of  the  first  parameter  and  the  process 
continued  until  it  is  found  that  the  temporary  best  values 
remain  unchanged  for  two  successive  cycles  and  can  be 
taken  as  the  best  values  for  the  exponential  parameters,, 
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INTEGRALS  FOR  THE  BASIS  FUNCTIONS  <*  AND 
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AIV-1  The  S  Matrices 


The  general  term,  ^ .  ,  for  the  matrix  is 

°  7  mm  7  — 


Smm '  ^2am  ^2%.  ^2am '  ^  1^2  > 

=  <%/^l/></^2^2 ft l/^2am£ 2/ ^  ^1  ^2am ' ?  2  “^l?2am‘^ 

“^11^22^  4im'^22  °12"Jmm '  '  ^mn1  nm'^ 

where  the  electrons  always  occur  in  the  order  1,2  and 
the  nuclei  as  1,2  and  S12  =<  V^Mpq.  =  ^ > 

^  -<$M  y  and  £.  .  is  the  Kronecker  delta0  Similarity 

1  J  *  f  J  10 


one  obtains 


4  2 

Dnn' 

,^i 

Drr ' 


S  2 

ss ' 


=  7?ll7?22[^nn'^ 

^11^22  “  S12V1]  £rr  ,7?  22 
^11^22  ”  °12V1‘1  ^ss '^11 


o  2 

11  °122mm '  °m&in '  %n '  °nm 


•  <Srm  1  ftr\m  '  ^ 


where  V±  =\2  =  slsl  +  &2&2' 


AIV-2  The  H  Matrices 

For  matrix  the  general  element  is 
Hmm'  =  <,^i^2am^2??l7?2/2ki:N  ^rkNSk*  V ^1^2am '^{^2^ 
-(^1  2^1 7  2//2k2N  ^"rkWSk  °  V  ^l£  2am  it  2/> 

=  “  H!m^2) 


These  two  terms  will  be  treated  separately  and  then 


.  . 

* 
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recombined.  For  the  first  term 

Hmm'  ^^m^l'^V  ^1AJA  +  ^1 

+  (  ^1^2am^2^1>?2//S2*  ^  %AZA  +^BIB"//^i^2am'^2>?l^^> 

=  ^2  <am/3/^t)[D^22<7^i/IA/^i>  +  db7?ii<^9?/2/IB//^2)-1 

+  4m'<£  2//S2/^2/>  *  [DAI?22<v?£  l^A^l^  +  DB^11 
<^'?2//IB//i^2/>  ^ 


where  A  and  B  refer  to  the  nuclei  and  1  and  2  to  the 
electrons  and 

°A  •<^U)/^.A/}i(1)>  ;  DA  =<^^(2)/4a/^2(2)'> 

DB  =^(1)/<T1b//^(1))>  ;  d2  =<^(2)/^a/^2(2)> 

■ ^  —Js  — ^ 

The  scalar  product  S^,  0  1^  can  be  written  as  ’  1^ 

=  StocInx  +  SkyIny  +  skzXkz  and  sinoe  the  bra-ket  nota" 
tion  has  put  the  functions  and  operators  together  correctly 

one  can  write 

,  (1)  =  H^,(lx)  +  H^,(ly)  +  ^,(lz)  ;  where  ,f or 
mm  mm  mm 

example, 

=>22  < VW>  [D1<VW?22 

+  DbWV?2>  ]  +  4m-<  ^2/Sx^2>K,f22<^l/Ix/7l> 
+  Db7h<'??2//Ix/72)  ] 

From  the  relations,  S  *  =  4|-  (?  ,  Sx  ^  =  4r°<,  and  defining 


mm 
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cT  ,  n  ,  one  obtains 


<iW0m  '>  =  T" 

<vw  =  2 


ID 


r,r 1 

1  h  h 

2  r  : 

„  2 

3  ,  s ' 

TT  V 

li2 

ln,n' 

T  ens: 

Inserting  these  integrals  gives 


■  T-^,»-^22dA  VSr ,r' hrhr '  4 ,r '+|lDBSs , s •  Vs ' 
*  V  +  C,m'  "2”  3n , n 1  ®=n®n 1  ^n,n'^?22^A  "5"  2r,r' 

Vr'^.r^ll^s.s'Vs^s.s-  ] 

1 

The  H'“;,(ly)  and  Hr  ,(lz)  contributions  may  be  obtained 
mm  mm 

in  a  similar  manner  and  combined  with  R'^.(Ix)  to  give 

mm  ° 


H^,(l)  =  i  ^  ,[dJ  ppWn  +  Di  Wp]  -  (-l)m'+1f  , 

mm  '  4  m,m  A  4  2  1  d  u  d  m ,  m 

*-DA  22w5  +  DB  11W6^  +  ^m,m'  22W1+  DB  llV 

+  V°A  22W3  *  DB  llV  +  U3(-IIA  22W5  +  DB  llV] 


where 

W  =  2  ,hh 
1  r ,  r 1  r  r 1  r , : 


W,  =  2^  L{  u.  .x  ,  «  , 

3  r,r  r  r  r,r  j 


^  h  ,  £ 


W0  =  2  ,  h  h  ,  f  , 

2  s ,  s  '  s  s  '  s  ,  s  ' 


(-l)s'+1  -p 

W„  =  21  t'  h  h  ,5 


W 


=  S^"1) 
5  r,r' 


r'+l 


hh , £  .  . 

r  r  r ,  r  ; 


4  s  ,  s  ' 

Wr  =  2^-1} 
6  s,s' 


s  s'  s , s 


s '  +1 


h  h  ,1 


s  s  '  s ,  s 
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U 


1  ~  2n ,  n '  ,  n ' 


!  U2 


s(-l)n’+1  £ 

L>-  -  I  bnfar,  I  n  n  « 


n,n' 


Jn^n  •  n ,  n 1 


TT  =  y(-l)n'  +1a.  f 

^3  ~  Sn,n'  -^n^n’  n,n'  ;  Uz 


=  2 


m ,  m 1  &m , m '  m,m' 


.  1 


U 


5  ”  ra,m' ' 


(-Dm'+V 


m,m 


U6 


^m ,  m ' ®m^m ' 


(-1) 


m '  +1 


4, 


m 


sin'  +1 

(  1}  S  ■ 

n 


^2  ”  ^nm 1  ^n^m 1  "^n  ,  m 1  ’  ^3  ^nm 1  1 

H^1,(2)  can  be  written  as 
mm 

1^,(2)  -^2VA,(yA 

+  ^  ^2AIA  *^2B1B  >,/  '^2^2am'>/l,?2y 

=  S12[  <am/S/^2^^2/c,m^  +  <'  om//f2^<-^2/'S/'0m^  ^ 
[of  ?22<?i/i7^>  ♦  ^2?n^2/f/^] 

where  ^2=<^/^/^ >  ;  D^2=^/i^/f|>. 

As  before  one  may  write 

H^,(2)  =  H^,(2x)  +  H^,(2y)  +  H^,(2z)  where,  for 
mm  mm 1  mm '  17  mm '  7 

example , 

=S12t<W?2 *W>  +<°m^2X;2/3x/om,>  ]* 
^A^VV*])  +  db2?ii<W*2>] 


which  becomes,  on  carrying  out  the  integrations, 

%  2 

^mm 1  ^ ^  “  “ZT  S12Sn,n' snsn*  ^^m,m' ^im'  +  ^nn'^n,m'] 
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•CDl2^f22Wl  +  3B2lllW2] 


12 


%  1 

The  H  ,(2y)  and  H'  ,(2z)  contributions  may  be  obtained 
mm '  °  mm 

in  a  similar  manner  to  give, 


■  £  Snot* 


4  °12L^n,n' ^n^n' ^  '  ^nm '  +  mn'  n,n'^ 


^12?22W1  +  DB2V?2)  -  Sn,n'Sn6n,(  +\,n'  ^ 

+  (-Dn'+1^n^n,m'>  ^22W5  +  Db\iV  +  2n, 

SnSn,(  (-i)a,+14nIrnjll.*(-i)",+1rIia4l.) 


(dA2^2W5  +  Db\iV] 

Combining  H^,(l)  and  3^7(2)  gives  finally ^ 

Hmm '  =  {^m,m'^22[I)A  ^22T,71  +  °B  ^  nW2]  +  (’1^m  1  ^m,m!^22 

[D^22w3  ♦  D^xW4]  +  (-!)“' +1^m,l22CBl^2»5+D^uW6] 

+  i'fUfD^22f/l  +  DI^1W2>  +  VDA?2^3  +  DB?llV 

+  VDA?22W5  +  DB^llV]}  “  I  S12i  Sn,n'snsn' 

+^n'C,m']  *  +  3b\iW2] 

-[(-Dn'+1^n,rnin,  ♦  (-i)"'+1^m, 

[d1%2w5  +  D^nw4]  +  C(-l)n,+14n,4B,+  (-l)m'+1 

CDA2?22W5  + 


But  it  can  be  readily  shown  that  W 


3’ 


and 


are  all 
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zero  and  H  *,  reduces  to 
mm 


Hram '  =  ^  ,  m  •  ^22CDA5?22W1  +  DB?ilW2]  +  (  ^  1^'m  m  • 

1*1  22W5  +  DB  llV  +  WUl[Df  22W1+DB  11W2] 

+  U3CDi  22w5+  DB  llV>  -S122n,n'SnSn- 
<-r&m,C,n'  +  [DA222W1+  DB211H2-1 

+  ^-Dn'+1  Sn-  L<  ^m’+1  in'  in'l 


db2?iiV)}  • 


In  an  analogous  manner  the  general  matrix  element  of  the 
other  matrices  maybe  reduced  to 

Jz  .  Jl 


im'  4 
2 


fe.a'  £l[Dl?2 2W1+  DB?11W2]  +  ("1)n+1  L>  11 

[DfW  DI?llV  +  £n'<VDI  22W1+  DB  11W2] 

+  VDfo>2W5  + 

^  ^nm'  ^n;n'  +^i ,m 1  4n' ]  l-DA2?<22Wl+  DB2;<11W2] 

+  [(-D“'+1  4m,  4n.  +a)m'+1  L. 

hiwe]  >} 

rr '  ‘  T-fe,r'^22CD1^22V  DjftiV  2Di2si2VlV2] 

*(-l)r’+1  fr,V D^22U6  +  D^liu3  -  2D12s12V1V3] 

+  4r-(VDB^22U4  +  D^11U1  “  2DfS12VlV  +  W6 


fDA  ^22S5+  DB 
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[Db^2U6  +  Db\iU3  2DB  S12V1V3^  J 

Hsf'  =  IR.S^II^bL^  +  Dl£lV  2DB2s12W 

+  (’1)S'+1  Is-  11  ^22U6  +  DB  ^11U3  -  2DB2s12T1V3] 

+  +  D!LU1  -  2D12s12V1V2] 

+  VD^22D6  +  DI\  1U3  -  2D!2si2VlV)J 


